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ABSTEACT 

This learning module is concerned with the 
temperature field, the heat transfer rates, and the coolant pressure 
drop in typical pressurized water reactor (PWE) fuel assemblies. As 
in all of the modules of this series, emphasis is placed on 
developing the theory and demonstrating its use with a simplified 
model. The heart of the module is the PWE Ther mal-Hydraulics Computer 
Code which solves for the radial temperature distributions in the 
fuel, cladding, and coolant of any axial station and then marches 
axially with an energy balance in the coolant. The code and its use 
are described in detail, including a listing and definition of all 
variables, a discussion of all input requirements and resulting 
output, an annotated flow-chart of the code, an explanation of all 
options in the code, and a listing of the code which gives enough 
comment statements to clearlv indicate the operational steps being 
performed. By proper specification of the options, the code can 
either be used as an individual entity to study thermal-hydraulic 
aspects exclusively, or as a subroutine in the total PWE module 
package to provide temperature feedback to the other modules. 
Examples are worked out using the code. (Author/SK) 
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THERMAL-HYDRAULICS MODULE > TH-1 

PRESSURIZED WATER REACTOR THERMAL-HYDRAULICS 

1.0 Object of Module 

The object of this module is to present the basic v [? of 
pressurized water reactor thermal-hydraulics . This rt:- \ 
demonstration of: 

• How the actual fuel geometry can be modeled for 
simplified thermal-hydraulic ainalysis. 

• What information is necessary to characterize the thermal- 
hydraulic behavior of the reactor. 

• The development of the theoretical relations that permit 
the computation of these thermal-hydraulic characteristics. 

• The actual calculation of this information for the reactor. 
This calculation requires the use of the PWR Thermal- 
Hydraulics Code, the description of which is included in 
this module. 

The thermal-hydraulic characteristics of the reactor are required 
for the determination of: 

• Fuel integrity 

• Cladding integrity 

• Coolant exit conditions 

• Pumping requirements 

9 Temperature feedback for reactor neutronics calculations. 
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2.0 Content of the Module 

This learning module contains the thermal-hydraulics of pressurized 
water reactors. Specifically, the module is concerned with the temperature 
field, the heat transfer rates and the coolant pressure drop in typical, 
pressurized water reactor fuel assemblies. 

As in all of the modules of this series, emphasis is placed on 
developing the theory and demonstrating its use with a simplified model. 
The model is carefully selected to insure that analyses based on it will 
exhibit all of the important thermal-hydraulic trends of the typical reac- 
tor. The description of the fuel for a typical pressurized water reactor 
and the modeling of the fuel which preserves the significant thermal- 
hydraulic characteristics are treated in the next section of this module. 

Following the geometry and modeling section, the basic theory 
governing the temperature distributions, heat transfer rates, pressure 
drops, and energy balance considerations is presented. The temperature 
distributions in the fuel and cladding are calculated assuming one-dimen- 
sional radial heat conduction. The pressure drop in the coolant channels 
and the heat transfer coefficient for use in Newton's law of cooling are 
calculated from empirical relations developed for reactor coolant channel 
flows. Energy balances for small axial segments of the coolant channel 
are used to step the solution in the axial direction. Simple examples,, 
illustrating the individual calculations, are worked out in detail for 
typical PWR conditions. 

The heart of the module is the PWR Thermal-Hydraulic*; Computer Code. 
Basically, the code solves for the radial temperature distributions in 
the fuel, cladding, and coolant at any axial station and then marches 
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axially with an energy balance in the coolant. The code and its use are 
described in detail. Included are a listing and definition of all vari- 
ables, a discussion of all input requirements and resulting output, an 
annotated flow chart of the code, an explanation of all options in the 
code, and a listing of the code which includes enough comment statements 
to clearly indicate the operational steps being performed. By proper 
specification of the options, the code can either be used as an individual 
entity to study thermal -hydraulic aspects exclusively, or as a subroutine 
in the total PWR module package to provide temperature feedback to the 
other modules. Examples are worked out using the code. In typical 
examples, the location and magnitude of the maximum fuel temperature in 
the PWR are found and the effect of oversized cladding on maximum fuel 
temperature and coolant outlet temperature are determined. 



EKLC 
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3.0 PWR Fuel Geometry and Its Model 

The typical pressurized water reactor core is approximately 
cylindrical in shape and is composed of 200 to 240 fuel assemblies. These 
fuel assemblies are rectangular in cross section, 12 to 14 feet long, and 
stand vertically in the core. Water gaps (the order of 0.20 in.) separate 
the fuel assemblies from one another. The fuel assemblies are contained 
inside a core shroud which in turn fits inside a cylindrical core support 
which is 10 to 12 feet in diameter. To provide a flattening of the neutron 
flux (and hence the power production) in the radial direction, the 
enrichment of the fuel is varied in cylindrical zones. 

The fuel assembly consists of 150 to 200 fuel rods, support grids 
holding the rods firmly at the ends, and spacers between the individual 
rods at several axial locations. All but about 1 foot of the 12 to 14 feet 
length of the fuel rods contain nuclear fuel and therefore generate heat 
internally during reactor operation. 

A cross-section of a typical fuel assembly is shown in Figure 1. 
The fuel rods are arranged in a square array within the fuel assembly. 
The coolant, which carries away the energy generateH within the fuel rods, 
flows in the spaces between the rods. Three types of coolant channels can 
be identified as shown in Figure 1. Thes^ correspond to the channels 
associated with central fuel rods, outer row fuel rods, and corner fuel 
rods . 

Noting that a great majority (144 out of 196 in a 14 by 14 array) of 
the coolant channels are central-type coolant channels, and that the outer 
row- and corner-type channels do not differ too markedly from the central- 
type (elongated about 10% in one or both directions, respectively) one of 



8 



- 6 - 



the central-type channels is selected as the model coolant channel. The 
model channel thus consists of a square with corners at the centers of 
the four fuel rods forming the channel minus the four quarter-circular 
segments of the fuel rods themselves. The ratio of center-to-center 
spacing, S, to the fuel rod diameter, D^, is typically 1.3. The hatched 
area in Figure 1 shows a model channel . 

Each flow channel serves segments of four fuel rods as can be seen 
in Figure 1. Since the four rods are identical and since the neutron 
flux does not vary significantly in the short distance from one rod to the 
next all four rods may be assumed to have the same operating characteris- 
tics. Thus the four fuel rod segments associated with each flow channel 
can be treated as one equivalent fuel rod. The power generation level in 
the model fuel rod will be taken as that of the average fuel rod in the 
core . 

The model fuel rod will have a variation in power density in the axial 
direction. This variation will be taken as 

q (Z) = q o cos ^ ^) 

unless a computed actual axial variation is provided from another module. 
? ? f 

In Equation 1, q represents the thermal source strength per unit volume 

? f t 

at any axial location Z, q ^ represents the thermal source strength per 

unit volume at the center of the fuel rod (Z=0) , and is the extrapolated 

1 1 1 fit 

height of the core. Both q and q are taken to be constant radially 

T I t 

throughout the fuel in the fuel rod. The magnitude of q ^ is representa- 
tive of that of the average fuel rod in the core. The extrapolated height 
is calculated from 

10 



H = H + 2 L (2) 
e e ^ ^ 

where H is the actual height of the reactor core and is an extrapolation 
length. The extrapolation length is the distance between the actual end 
of the core and the location where an extrapolation of the waveform re- 
presenting the actual neutron flux distribution within the core goes to 
zero. From neutron diffusion theory can be shown to be about one 
migration length for a bare core. The migration length can be calculated 
from neutron diffusion theory [1,2,3]. 

The PWR fuel rods are 12 to 14 feet long and about 0.45 inch in 
diameter. About 8 to 10 inches of the total length is used for end 
fittings and a fission gas plenum. The remaining length contains nuclear 
fuel, usually uranium dioxide. The fuel is typically in the form of cylin- 
drical pellets with a diameter of about 0.39 inch. These fuel pellets 
are stacked inside a tube of cladding, usually a zirconium alloy with 
about a 0.025 inch thick wall. The 0.003 to 0.005 inch gap between the 
fuel and the inside of the cladding is pressurized with high-purity 
helium. 

For simplicity, the helium gap is neglected in the model of the fuel 

rod. The model fuel rod thus consists of fuel of diameter contained in 

a tube of cladding with an inside diameter D- and an outside diameter D . 

r s 

For the model fuel rod, the length of the end fittings will each be 
assumed to be 1% of the active fuel length and the length of the gas 
plenum on the top of the fuel will be taken as 5% of the active fuel length, 



*Numbers in brackets refer to items in References. 
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4,0 PWR Thermal-Hydraulic Tlieory 



A.l Internal Heat Generation 

As a result of nuclear fission in the fuel, heat Is generated. The 

rate of energy generation in the fuel per unit volume is called the 

I I I 

"volumetric thermal source strength," q , and can be calculated from 

q'" - E^^^ 0^ (() (3) 

where is the energy released per fission reaction (energy dimensions), 

N^^ is the fissionable fuel density (fissionable nuclei per unit volume), 

is the effective fission microscopic cross section (dimensions of area), 

and ({) is the neutron flux (neutrons per unit area per unit time) . Note 
I 1 I 

that q has dimensions of energy per unit volume. The details of this 
calculation of volumetric thermal source strength are found in Reactor 
Statics Module 8. 

The volumetric thermal source strength varies throughout the reactor 

since ({), and perhaps N^^, vary. For a cylindrical reactor, axial symmetry 

is generally a reasonable approximation. Therefore q reduces to a 

function of only the radius, R, and the axial position, 2. Across any 

I I I 

single fuel rod, the change in q is small (due to small change in R) 
I 1 1 

and q can therefore be considered constant across Its cross section. 

I f f 

Thus for any single fuel rod q is a function of only the axial position. 

Also for any small axial segment of a fuel rod (the order of 1% of the 

f f f 

total active length) the change in q is moderate. Using a constant 
average value for any such segment thus introduces little error into the 
analysis. This approximation is used in this module. 
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The heat generation in any small segment of the fuel can be obtained 
by multiplying the volumetric thermal source strength by the volume of the 
segment. For steady state conditions this energy must be removed from 
the fuel. The mechanism by which this heat transfer occurs within the 
fuel is thermal conduction. 
Example 1 

A PWR containing 25,000 fuel rods generates 2000 flw(t). The cylindrical 
core of the reactor is 10 ft long and 7 ft in diameter. The fuel consists 
of 0.3 in. diameter fuel pellets contained in cladding which has 
an outside diameter of 0.35 in. The fuel rods are spaced 0.45 in. center- 
to-center in a square array. The volumetric thermal source strength of 

1 1 1 I I 1 

the fuel varies as q ^ q cos(nZ/H^). The extrapolation length for 

I 1 1 

the fuel is 6 in. Find the power per unit length (in kw/ft) and q ^ for 

the average fuel rod in the core. 

Solution 

The power generated in the average fuel rod is 

P =P_I£M1 . 2000 10^ 80,000 watt = 80 kw. 
ave N /3,UUU 

, , . P 80 kw _ Q kw 

Thus, per unit length the power is — = ToTit ~ "ft* 

The power produced by the average fuel rod can also be calculated by 
integrating the power produced in each differential volume of the fuel 
rod. Thus, 



P 
ave 



1 1 1 

q dV « 



rH/2 

III 2 



V ^H/2 



q cos (nZ/H )(nD"/4)dZ 
o e r 



H/2 



= (H ul q"'/2) sin(nU/2H ). 
e r o e 
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t I I 

Solving for q 



I 1 1 



2P 

ave 

° II sin(nH/2H ) 
e f e 



(2) (80.000 watt) (3.413 ;;^^^^) 



(11 ft) (0.3 in.)^ i^^.^ sin[(n)(10 ft)/(2)(ll ft)] 



= 8.025 X 10^ Btu/hr-ft^. 
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4.2 Temperature Distribution in the Fuel 

For conduction heat transfer the heat flux Is proportional to the 
normal temperature gradient, \7hen the proportionality constant is inserted 

^n n dX^ 

where is the heat transfer rate in the n-direction (energy / time) , 3T/5X^ 
is the temperature gradient in the n-direction (degrees/length) , is the 
area normal to the n-direction, and k is the thermal conductivity of the 
material (energy /degree-leng th-time) . Equation 4, which is Fourier's law 
of heat conduction, relates the heat transfer to the temperature field 
and is also the defining equation for the thermal conductivity. The ther- 
mal conductivity is a material property and its magnitude in general varies 
with the temperature of the material. Heat transfer properties of various 
reactor materials are tabulated as functions of temperature in References 
4, 5, and 6. To reduce the complexity of heat transfer calculations, the 
thermal conductivity is often assumed to be constant and evaluated at an 
average temperature. The minus sign in Equation 4 assures that the 
heat transfer is in the direction of decreasing temperature. Equation 
4 shows that temperature gradients are required for heat transfer^ In 
nuclear reactor applications, where there are high heat transfer .-ates, 
large temperature variations occur. One of the primary tasks of reactor 
thermal-hydraulic analysis is the prediction of this temperature distri- 
bution in the fuel. 

The heat conduction equation in cylindrical coordinates 

i!l + ill + 1.1^ + 1^ + slLL = i 11 

3r2 ^ r2 302 32^ ^ ^ ^ ^ 
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and the initial and boundary conditions prescribe the temperature 
distribution within the fuel rod. In this relation, the thermal conductivi- 
ty, k, and the thermal diffusivity, a, have been assumed constant. The 
development of this equation is found in heat transfer texts [4, 5, 6] 
and several simplified cases are given as exercises for the student. 
For steady sta^te "conditions , the unsteady term on the right hand side of 
the equation is zero, and for axial symmetry, the 0 variation disappears. 
It may also be observed that since the length of a fuel pin is much greater 
than its radius, the temperature gradients in the radial direction will be 
much greater than the temperature gradients in the axial direction. 
Therefore, to a good approximation, the heat transfer in the axial 
direction can be neglected with respect to that in the radial direction, 

and the resulting^, temperature distribution and heat transfer reduced 

1 1 1 

to a one-dimensional case for any axial segment in which q may be 
assumed constant. The governing differential equation for this case 
reduces to the ordinary differential equation 

.dr^ r dr 

The solution of Equation 6 subject to the boundary conditions 

T = T at r = 0, 

(7) 

4^ = 0 at r « 0 

yields the temperature distribution in the fuel. The second boundary 
condition is obtained from the observation that the temperature distri- 
bution must be continuous across the center of the cylinder. Observing 
that both q and k are constant and that 
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dr2 ^ r dr I dr (8) 



Equation 6 can be written as 



d_ 
dr 



r P U - ^ 



I I I 



dr I (9) 



Integrating twice results in 



T = - •? 7- + C, In r + C,. 

^ ^ 2 . (10) 

Applying the boundary conditions of Equation 7, the integration constants 
are 



= 0, 



C, = T . 
2 o 



(11) 



Substituting into Equation 10 gives 

^ o 4kf ^ • (12) 

This relation shows the temperature distribution in the fuel to be 
parabolic with maximum temperature at the center- The heat transfer rate 
through any cylindrical shell can be calculated from Fourier's law which 
takes the form 

r f r dr 

where = 2'iTr(AL), AL being the length of the cylindrical shell. 

Of particular interest are the temperature and heat transfer at the 
surface of the fuel; i.e., at r = r^ = At this location 

f f f 



f O 4kp ''f ' (14) 
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dT 



r = 



irrf^(AL) q 



1 1 1 



(15) 



Noting that irr^ AL is the volume of the fuel rod segment of length AL, it 
is observed that the heat transfer out of the r = r^ cylindrical shell is 

I I f 

indeed equal to the total energy generated as calculated from (q ) 
(fuel volume). For one-dimensional heat transfer all of the energy 
generated within the fuel must be transferred out through the surface. 
Example 2 

Equation 14 relates the temperatures at the center and surface of a 
heat-generating cylindrical fuel rod. This equation contains the thermal 
conductivity of the fuel which, in general, is a function of temperature. 
However, the analysis leading to Equation 14 assumes the thermal conduc- 
tivity to be constant. Determine the temperature at which to evaluate 
to make the assumption of constant k^ compatible with a k^ that varies 
linearly with temperature. 
Solution 

The method of solution is to develop an expression relating T^ and 

T which assumes k^ - a + bT and then compare this result with Equation 
s f 

14. We begin with Fourier's law of heat conduction, 

dT 

^ f rdr 

I 1 1 



Recognizing that q = q (volume) and - 27rr(AL): 



I I I 9 dT 

q 7rr (AL) - -k^ 27rr(AL)^. 
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Simplifying, separating variables, and substituting for k^; 



q rdr = -2(a + bT)dT . 

» I \ 

Integrating from r = 0 to r = r^ with q , a, b constant, 

;nr^ j— — iT, 



2 



= -2 



aT + I 



2 

1 1 I h 2 2 



= [a + I (T^ + T^)] (T^ - Tj) 



III 2 



T - T- = 
o f 



4[a + |(T^ + T^)] 



Comparing this with Equation 14 written as 



111 2 



T - = -7-;- 
o f 4 k. 



we see that the relations are equivalent for 



k^ = a 4- - (T^ 4- T^) 



Note that this is precisely equal to k evaluated at the arithmetic 

T + 

average fuel temperature; i.e., k( ^ )• 



ERIC 
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4, 3 Temperature Distribution in the Cladding 

The temperature distribution in the cladding can be calculated by 

methods similar to those used for the fuel. Noting that no energy is 

generated in the fuel, the governing equation simplifies to 

f (r 4^) = 0. (16) 
dr dr 

The boundary conditions obtained by matching temperature and heat transfer 

at the fuel-cladding interface are 

T = T^ at r =» r^, 

ii= qf ^fq'" ^ (17) 

dr " 2TTk r.(AL) "2k ^ ^f 

c t c 

where k is the thermal conductivity of the cladding, 
c 

Integrating Equation 16 twice yields 

T « C3 In r + C^. (18) 



Applying the boundary conditions in Equation 17, 

rW (19) 

r = T^ + — 1^ 1^ 

^4 ^f 2k^ ^ 

Substituting into Equation 18 yields the temperature distribution in the 

cladding, 

f 1 r 



T = T. - -^r In . (20) 



Or, if q is eliminated in favor of q^, 

*rhe temperature at the surface of the cladding, T^, is obtained by 
evaluating Equation 21 at r^. 
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f t s 



^ " 27Tk^(AL) • (22) 

Note also that the same heat transfer rate must occur through each 
cylindrical layer for steady state conditions; i.e., Qg* 
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4.4 Temperature Drop in the Convective Layer 

All of the heat generated within the fuel ultimately must be 
transferred to the coolant. Thus all of this energy must be transferred 
through the fluid layer near the outer surface of the cladding. The 
heat transfer mechanism, wherein the energy is carried away or convected 
from the solid surface by a fluid in motion, is called convection. The 
heat transfer rate for convective heat transfer is related to the tempera- 
ture difference between the surface and the bulk fluid, the driving 
potential for the heat transfer, by Newton's law of cooling, 

qs"^^s ('^s- V- (23) 

This relation may also be taken as the defining equation for the heat 
transfer coefficient, h. The bulk temperature T^, is a mass-weighted 
average temperature of the fluid in the flow channel. It is formally 
defined by 

fpCUT^ 

This is the temperature that a thermometer would indicate if iimersed in 
a cup of fluid collected from the discharge of the flow channel at the 
given location. 

Relations for calculating h from the flow characteristics and 
coolant properties have been developed. These will be treated in a later 
section of the mocule. For the moment, the quantity of interest is the 
temperature drop across the convective layer. In terms of h, this can be 
obtained from Equation 23 as 
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" ""^B " 2irr^(AL)h (25) 

where is the outer radius of the cladding. 
Example 3. 

Consider the PWR described in Example 1. The bulk temperature of 
the coolant changes from 491 F to 609 F as the coolant flows through the 
reactor. Given that the heat transfer coefficient is 6000 Btu/hr-f t^-F, 
the thermal conductivity of the fuel and cladding are 2.40 and 6.90 
Btu/hr-ft-F, respectively, determine the fuel centerline temperature at 
the axial center of the core. 
Solution . 

Half of the total energy transfer to the coolant occurs in the first • 
half of the core. Thus the coolant bulk temperature at the axial center 
of the core is 550 F. To this, we must add the temperature rises across 
the convective layer, the cladding, and the fuel to obtain the fuel center- 
line temperature. 

The heat transfer rate through the convective layer at 2 = 0 is 

/ III o 

"^s ~ o ^^"^-^ volume) = q ^ irD^(AL)/4 

= (8.025 X 10^ Btu/hr-ft^)(n/4)(0.3 in.)2(ft/i2 in.)^(LL ft) 
= 39393 (AL) Btu/hr . 

From Equation 23, 



T - T = — = 39393 AL Btu/hr 

s (6000 Btu/hr-ft^-F)n(0.35 in.) (ft/12 in.)(AL ft) 

- 71.6 F. 
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The heat transfer rate through the cladding is the same as that 
through the convective layer. Equation 22 yields 

_s_ 

" ^s " 2TTk (AL) r. 

c t 

^ (39393 AL Btu/hr) ln(0. 175/0. 15) 
2tt(6.90 Btu/hr-ft-F) (AL ft) 

= 140.1 F. 

From equation 14 the temperature rise through the fuel is found as 

^ ^ _ (8.025 X 10^ Btu/hr-ft^)(0.15 in.) (ft/12 in.) 

"^o " ^f ' ~4~k^ 4(2.40 Btu/hr-ft-F) 

= 1306.2 F. 

Adding these temperature rises to the bulk temperature, 

= ^^o - ^f^ ^ ^^f - "^s^ + ^^s - h'> + ^B 
= 1306.2 + 140.1 + 71.6 + 550 - 2067.9 F. 
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4.5 Pressure Drop in the Coolant 

All of the energy generated by nuclear fission in the fuel must be 
carried out of the reactor by the coolant which flows in the channels 
between the fuel rods. In a PWR the coolant is water which is force- 
circulated by a pump. Enough pressure head must be provided by the pump 
to overcome the pressure losses in the water flow loop. The entire loop 
consists of the reactor flow channels, the heat exchangers, water quality 
maintenance components and connecting piping. Of interest here is the 
determination of the pressure loss incurred in the flow through the 
reactor core. Two types of pressure losses will be considered. These are 
the frictional pressure loss along the coolant channel and the entrance 
and exit permanent pressure losses. 

The frictional pressure loss is a manifestation of the shear stress 
on the flowing fluid by the walls of the flow channel. This pressure loss 
is calculated from 

AP = f ^ ^ 
^ • (26) 

This relation may also be interpreted as the defining equation for the 
friction factor, f. The bulk velocity, U^, found in Equation 26 is the 

D 

velocity averaged across the flow channel. It is defined by 



"b = I / U dA. (27) 
A 

This velocity is related to the mass flow rate by 

m = pAUg (28) 
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or to the mass velocity by 



"h- (29) 



The equivalent diameter, D^, is used to relate experimental data from 
circular tubes, which is abundant in the literature, to different geometries 
The equivalent diameter is calculated from 

4A 

I>e " • (30) 

2 

Note that for a circular cross section the flow area is IID /4 and p = ITD, 

w 

Thus for a circle the equivalent diameter is equal to the actual diameter, 
as required. This> of course, is why the definition of Equation 30 is xised 

The value of the friction factor depends on the flow acid surface 
conditions in the channel. The flow conditions are characterized by the 
Reynolds number, 

" "T" • (31) 

It is well known that for Re below a critical value the flow is laminar, 
and for Re above the critical value the flow is turbulent. For internal 
flows this critical Re is about 2000. Laminar flow may be thought of as 
an ordered process in which fluid layers slide over one another, being 
retarded only by the molecular interaction between the layers. The vis- 
cosity of the fluid quantifies the magnitude of this interaction. In 
laminar flow any disturbance in the fluid is damped by the viscous action. 
In turbulent flow there is an additional random transport mechanism 
operable. This mechanism may be modeled as eddies (finite sized patches 
of flxiid which retain their characteristics for finite times) moving 
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throughout the fluid, transporting mass, Tnomentum, and energy by virtue 
of their movement. This action is quite violent and results in transport 
rates much greater than those by purely molecular activity in laminar flow. 
In turbulent flow, disturbances in the flow field grow and propagate 
resulting in increased turbulence levels downstream of the disturbances. 
Even in highly turbulent flows, there exists a layer near the solid bounding 
surfaces where the presence of the wall retards the penetration of eddies 
and therefore acts as a flow stabilizer. This results in a laminar sub- 
layer existing near any solid boundary. Even though this layer is very 
thin, it is very iii5)ortant, since much of the temperature and velocity 
change between wall and bulk conditions occurs in this layer. In reactor 
coolant channels the flow is highly turbulent with Re of the order of 
100,000 and greater being quite common. 

The surface condition of the flow channel is characterized by the 
ratio of the surface roughness height to the equivalent diameter of the 
flow channel, e/D^. The flow channel exhibits smooth tube behavior if the 
roughness height is less than the thickness of the laminar sublayer. For 
reactor coolant flow channels the surface conditions are controlled so 
that this criterion is met. 

The calculation of the friction factor for the reactor coolant flow 

channels is a two-step procedure. First, the smooth circular tube friction 

factor is determined, and then a correction factor is applied to account 

for the actual channel shape. The smooth circular tube friction factor 

can be determined from a Moody chart where f . is plotted versus Re. 

clrc 

These charts are found in most fluid mechanics texts and handbooks [7, 8, 
9, 10]. An alternate method is to calculate the friction factor from a 
correlation in equation form. One of the most widely accepted for reactor 
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work Is 



f = 0>184 (32) 

^circ ^ 0.2 • ^ ^ 

Re 



This relation, like the Moody diagram, was obtained from a curve fit of 

experimental data in smooth circular tubes. Tlie friction factor defined 

in this module is the Darcy-Weisbach friction factor. Care must be taken 

not to confuse this fraction factor with the Fanning friction factor (fp^^ 

f /4) also found in the literature. A correction of the circular tube 
D-W 

data to the rod bundle geometry has been obtained by Deissler and Taylor 
[11]- Their results, expressing f/f^^^rc ^ function of S/D^, can be 
approximated by the following polynomial fit: 

s s 2 s 3 

f/f . =-3.475 + 8.053 ^ - 4.705 ) + 0.9162 (tt-) . (33) 
circ Dg Dg Dg ^^^^ 

The permanent pressure losses at the inlet and exit are the result 
of increased viscous energy dissipation resulting from increased turbulent 
activity. These losses may be calculated from 

In this relation, K is the resistance coefficient for the inlet or exit 
geometry. The inlet may be approximated as a sudden contraction from a 
very large diameter to the equivalent diameter of the flow channel. For 
this sudden contraction, the resistance coefficient is 0.5 [10]. 
Similarly the exit may be approximated as a sudden expansion from to 
a very large diameter. For this expansion, K is 1.0 [10]. 
Example 4. 

The PWR described in Example 1 has an inlet coolant temperature of 
491 F, an inlet coolant pressure of 2250 psia, and a coolant mass velocity 
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6 2 

of 2.5 X 10 Ibin/hr-ft • The coolant exit temperature is 609 F, 
Estimate the frictional pressure loss in the core. Neglect entrance 
and exit losses. 
Solution . 

The coolant properties are evaluated at the average coolant temp 
ature. At 550 F, 2250 psia: 

U = 0.2336 Ibra/hr-ft 
p = 46.47 Ibin/ft^ 
The average coolant velocity is found from 

= G/p = (2.5 X 10^ lbm/hr-ft^)/(46.47 Ibm/ft^) = 53,800 ft/hr. 

The equivalent diameter is computed next. 

A = S - y I) 
r 4 s 

= (0.45 in.)^ - -J (0.35 in.)^ = 0.1063 in.^ 

p = IID = 11(0.35 in.) = 1.100 in. 
"^w s 

= 4 A^/p^ = 4(0.1063 in.^)/(1.100 in.) 
= 0.386 in. = 0.0322 ft. 
The Reynolds number is 
Re - p/p 

= (53,800 ft/hr) (0.0322 ft) (46.47 Ibm/f t^)/(0.2336 Ibm/hr-ft) 
= 344,600. 



Using a circular tube friction factor correlation, 

f , « 0.184/Re°-^ - 0.184/(344,600)°-^ 
circ 

= 0.0144. 



Applflring the correction for the actual reactor flow channel, 

/ • 

f = f [- 3.475 + 8.053(S/D ) - 4.705(S/D^)^ + 0.916(S/D )^] 

circ s s o 

2 3 

= (0.0144) [- 3.475 + 8.053 (§^) - 4.705 (§^) + 0.916(|y||) ] 



= 0.0151. 

From the defining equation for the friction factor, 
APy = f(AL/D^) p U^/2g^ 

= (0.0151) (120 in. /O. 386 in.) (46.47 Ibm/ft^) (53,800 ft/hr)^ 

2 2 
-/(3600 sec/hr) (2) (32.2 Ibm-ft/sec - Ibf) 

= 756.5 Ibf/ft^ 

= 5.25 psi. 
Example 5^^ , 

Each coolant ch'^annel in a fuel assembly has the same pressure drop. 
If all coolant channels are not identical the flow will redistribute such 
that this equal pressure drop is attained. Determine the relation 
betw^n average velocity in the coolant channel (or equivalently the 



maS^velocity)^ and the .equivalent diameter of the channel which governs 



the flow distribution. 
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Solution. 



The frictional pressure drop is given by 

,2 



AP = f 



Substituting for the friction factor. 



AP = 



0.184 AL ^ 
e 



0.184 AL 



0.2 D 2g 

e c 



D e 



if constant density and viscosity are assumed. For channels 1 and 2, 
each having equal pressure drop. 



or 



1.8, 1.2 ^ 1.8, 1.2 
^Bl 'el B2 ' e2 



2/3 
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A. 6 The Heat Transfer Coefficient 

The heat transfer coefficient for reactor coolant channels is 
determined by a two-step procedure similar to that used previously for 
the friction factor. The heat transfer coefficient for flow in circular 
tubes is found first and then a correction is applied to obtain h for the 
actual reactor coolant channel geometry. 

No direct analytical methods for predicting h for turbulent flow in 
circular tubes exist. However, by the use of an analogy between heat and 
momentum transport, the dependence of h on the flow field characteristics, 
i.e., on the Re can be predicted. Then, by applying the results of a heat 
transfer analysis for laminar flow over a flat plate, an tstimate of the 
Prandtl number dependence on heat transfer can be obtained. The 
Pr may be interpreted as a dimensionless modulus relating the temperature 
field in the fluid to the flow field. For Pr of unity, similar velocity 
and temperature profiles exist in the fluid. Finally, the predicted 
expression for h is compared to empirical correlations found to fit the 
available experimental data. 

Tlie heat transfer across a fluid layer in laminar flow may be 
calculated from Fourier's law, 

dT 

q - - k A ^, (35) 

where y is the direction normal to the fluid flow. Similarly, the shear 
stress in the fluid can be related to the velocity gradient by 

X = . g (36) 

which defines the viscosity coefficient y. These relations can be 
rearranged to yield 

q dT 
PCpA °y (37) 
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and 



T _ .. dU 

(38) 



P dy 



k 

In these relations, a - --r- is the iiK>lecular diffusivity of heat and v 
is the kinematic viscosity or the molecular diffusivity of momentum. We 
now postulate that the eddy transport of heat and momentum present in 
turbulent flow can be expressed in the same form. To this end, we define 
the eddy diffusivity of heat, e^^, and the eddy diffusivity of momentum, e^, 
as the parameters, that, when multiplied by the appropriate gradient, 
yield the corresponding transport rate for turbulent flow. In general e^^ 
and vary throughout the flow field. For the combined molecular trans- 
port plus turbulent eddy transport, 

^ = - (a 4- sj ^ 



^CpA " ^" ^ V d^F (39) 
and 

5 = + ^M^ d^- 

(40) 

We now assume that heat and momentum are transported by analogous 

processes and at the same rate. This requires that a = v and e_, = e^. 

H M 

Since the Pr = v/a, the first condition is equivalent to Pr = 1. Similarly 

a turbulent Prandtl number may be defined as Pr = e^7£„. Therefore the 

t M H 

second condition is equivalent to Pr^ = 1. The basic assumption we have 
made also implies that both q and t vary in the same way across the field; 
i.e. , 



C„At 



Constant 



p"' SVw' (41) 
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if the constant is expressed in terms of quantities at the surface or wall. 

Dividing Equation 39 by Equation 40 and using the assumption dis- 
cussed above. 



dU = -dT. 



C A T„ (42) 
p w w 



Integrating this expression between wall conditions and fluid bulk 
conditions 9 ^ '^B 



dT. 



CpVw (A3) 
w w 



Since U = 0, this reduces to 
w 



w"b „ _ t (aa) 



Cp^w^w 



To introduce h into this expression, Newton's law of cooling (Equation 23) 
is used to give 



% 



Substituting Equation 45 into Equation 44 , 



. (46) 



\ ' 

The shear stress at the wall can be eliminated from this expression 
in favor of the friction factor. To do this, is first related to the 
pressure loss over the length of the pipe, AL, and then Equation 26 re- 
lating the pressure loss and the friction factor is used. Consider the 
fluid in a section of pipe of diameter, D, and length, AL. For fully 
developed flow, the sum of the forces acting in the flow direction must 
equal zero. The forces that must be considered are the pressure force on 
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the upstream circular face, |- D P, the pressure force on the downstream 
face, (P-AP) , and the shear force along the periphery , AttD Lt^. Summing 

these forces with proper regard as to their direction gives 



y D^P = 7T DALt + t ^ (P-AP) 

and 

w 4AL 

Using Equation 26 for AP, this becomes 

-_D_ f AL ^ 
""^w 4AL D 2 

Substituting Equation 48 into Equation 47, 



(47) 



(48) 



(49) 



Introducing the Stanton Number, a convenient dimensionless grouping common 
in heat transfer work. Equation 49 becomes 

St = — \ = 4 . (50) 

pCp Ug 8 

Using Equation 32 for f, 

-0.2 

St = 0.023 R„ 

(51) 

This result is often found expressed in terms of the Nusselt number, defined 
as 

(52) 

One common interpretation of the Nu is that it is the ratio of the actual 
convection heat transfer from a surface to the heat transfer assuming that 
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only molecular conduction in a stagnant fluid is present. Noting that 

h . ^ - Nu 

" ^ ^ ^ CpU PDI^ pr Re 

^ k ^ ^ y 

we can rewrite Equation 51 as 

Nu - 0.023 Pr Re°-^. ^^^^ 

The Pr dependence found in Equation 54 is not correct. One 
reason why this incorrect dependence appears is that Pr » 1 was assumed 
early in the development; therefore one would not expect the Pr dependence 
to be properly represented in the result. The analysis of laminar flow 

heat transfer from a flat plate (one of the few cases which can be treated 

1/ 3 

analytically) yields a Pr dependence in an expression for the Nxisselt 

1/3 

number. In practice it turns out that this Pr dependence is reasonably 
accurate for the turbulent pipe flow case of interest here as well. In 
fact some of the empirical correlations found in the literature exhibit 
this Pr"^^^ dependence in the Nu(Pr, Re) relation. However, the relation 
that is the best known, most widely used, and is recommended here is the 
Dittus-Boelter Equation [12], 



0.8 „ n 

Nu - 0.023 Re" Pr • 



(55) 



For the fluid in the tube being heated, the recoimnended value of n is 0,4 
and for the fluid in the tube being cooled, the recommended value of n is 
0,3. The former case is the one found in reactor work. In Equation 55 
several temperature dependent fluid properties appear (k, p, y, C^) . 
These properties are to be evaluated at the average bulk temperature of 
the fluid for the segment of tubing of interest. 



36 



- 33 - 



The heat transfer coefficient for circular tubes calculated from 
Equation 55 is then related to the heat transfer coefficient for the 
reactor coolant channel. This is done by using the equivalent diameter 
of the coolant channel in the Nusselt number and Reynolds number in the 
circular tube Nusselt number correlation and then applying the correc- 
tion of Deissler and Taylor [11]. The Deissler and Taylor correction is 
the same as that used for the friction factor; i,e, it may be approximated 
by 

N 23 

- -3.475 + 8.053 I- - 4 . 705 (|-) + 0.9162 (^) . (3,, 

The heat transfer coefficient determined from this relation is the average 
value around the periphery of the fuel rod. 
Example 6 . 

Determine the average heat transfer coefficient for the PWR described 
in Examples 1 and 4. 
Solution . 

The average heat transfer coefficient corresponds closely to that 

calculated at the average coolant conditions; i.e., 550 F and 2250 psia. 

Additional fluid properties required for this calculation are 

C = 1.258 Btu/lbm-F, 
P 

k = 0.3346 Btu/hr-ft-F. 

The Prandtl nxiniber is 

Pr = C u/k 
P 

= (1.258 Btu/lbm-F) (0.2336 lbm/hr-ft)/0.3346 Btu/hr-ft-F) 
« 0.878. 

From the circular tube Nusselt number correlation, 

Nu , = 0.023 Re°'^ Pr°'^ 
circ 

= 0.023 (344,600)°*^ (0.878)°'^ = 587.5. 
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Applying the correction for the actual flow channel, 



Nu = Nu , [-3.475 + 8.053(S/D ) -4.705(S/D ) + 0.916(S/D ) ] 
circ s s s 



E 



= (587.5) -3.475 + 8.053 



,0^ 
0.35 



-4.705 



0.45 



0.35 



+ 0.916 



0.45 



0.35 



615.7 . 



Calculation of h from the definition of Nu gives 



h = (Nu) k/D 



= (615.7) (0.3346 Btu/hr-ft-F)/ (0.0322 ft) 



= 6397 Btu/hr-ft-F . 
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4.7 Incremental Energy Balances 

As heat is transferred from the fuel rod to the coolant, the 

temperature of the coolant rises. Thus, the coolant temperature varies 

from a minimum at the inlet of the coolant channel to a maximum at the 

exit of the fuel bundle. The bulk temperature at any axial position can 

be determined from an energy balance where the total energy added by heat 

transfer to the coolant is equated to the energy rise of the coolant, i.e., 

Z 



mi = m±^^ J % (Z) dz. (57) 
-H/2 

Recognizing that for subcooled water Ai = C^AT^ and solving for the local 
bulk temperature, this becomes 

Z 

■ . ^ (Z) dZ. 

Equation 58 can be used to determine the bulk temperature increase for any 
segment, AL, of the coolant channel by 



Z +AL/2 
c 



(59) 



where is the center of the increment under consideration. Defining an 

_ t I _ 

average surface heat flux for the increment, q = q /A , the incremental 

s s s 

bulk temperature rise becomes 

TTDgAL _,, 

^B2 " ^Bl = ^b-C~ % • (60) 
P 

II 

The local surface heat flux, q^CZ), is related to the local volumetric 
thermal source strength by observing that all of the energy generated in the 
fuel is transferred out through the surface of the cladding; i.e.. 
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TTDg dZ qg'(Z) = J d2 dZ q"'(Z). ^^^^ 



This simplifies to 



Similarly, 



L '"r7.>>. (62) 



. ^ q . (63) 
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4.8 Assessment of Local Coolant Conditions 

Two important checks on coolant conditions should be made to insure 
that the reactor is adequately cooled. The first is a determination of 
whether the fluid bulk is indeed subcooled. As the coolant flows through 
the reactor, its energy is increased by heat transfer and its pressure is 
decreased by flow losses. Both of these effects reduce the subcooling of 
the fluid. To insure bulk subcooling, the local bulk temperati7re must be 
less than the saturation temperature at the local pressure. The critical 
location, of course, is at the exit of the channel, where the coolant 
bulk temperature is the maximum and the pressure is at a minimum. 

The other determination that should be made is whether there is any 
local boiling near the surface of the fuel rod. The criterion for this 
determination is whether the maximum coolant temperature at any cross 
section is greater than the local saturation temperature. The maximum 
coolant temperature at any cross ? . "ion occurs in the fluid layer immed- 
iately adjacent to the cladding is magnitude approaches the surface 
temperature of the cladding. The axial location of the maximum cladding 
surface temperature can be obtained by writing an expression for T (Z) 

HI 

in terms of q ^ , H^, D^, D^, m, C^ and h, differentiating with respect 

to Z, setting the derivative equal to zero, and solving for Z . The 

m 

maximum cladding surface temperature can then be obtained l.y evaluating 

T (Z) at Z . The details of this calculation are left as an exercise for 
s m 

the student. 
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4.9 Subcooled Flow Boiling 

Even though the bulk fluid is subcooled for PWR operation the 
possibility of local or subcooled flow boiling exists. Whenever the outer 
cladding temperature is greater than the local saturation temperature sub- 
cooled boiling may occur. Such conditions are obtained in PWR'S during 
peak loads. Under these conditions there is a possibility of boiling 

burnout of the surface. To assess this possibility, the actual heat flux 

II " 
at the surface, q^, is compared to the critical heat flux, q^, computed 

for subcooled flow boiling. The critical heat flux (CHF) , also, called the 

burnout heat flux or departure from nucleate boiling (DNB) , corresponds to 

the heat flux at which blanketing of the heat transfer surface by coalesced 

Vapor bubbles is just inuninent. Any further increase in heat flux would 

require the heat transfer to occur through the vapor blanket. Because of 

the low conductivity of the vapor blanket, a catastrophic increase in 

surface temperature would result and generally lead to melting of the 

cladding. 

One of the most widely used subcooled flow boiling critical heat 
flux correlations is the ANL correlation of Jens and Lottes [13]. This 
correlation is based on a large number of experimental data for subcooled 
water flowing vertically upward in a flow channel under a range of 
operating conditions. The conditions spanned by the data are mass velocity 
from 0.96 x 10^ to 7.8 x 10^ Ibm/hr-ft, pressure from 500 to 3000 psia, and 
subcooling from 5 to 163 F. The correlation for the critical heat flux, 

" . c (-^)'" (T - TJ°-". (64) 
c ^ \q6^ ^ sat ^B^ 

contains two pressure dependent "constants" which are given in Table 1 
below. 
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TABLE 1 

The "constants" C and m for Jens and Lottes Correlation 

W'ater Pressure, 

psia C m 



500 0, 817x10^ 0,160 

1000 0 ,626x102 0,275 

2000 0.445x10^ 0.500 

3000 0.250x10^ 0.730 



Note that the Jens an^ Eottes correlation breaks down as the 
subcooling approaches zero. For this case the more complex Westinghouse 
correlation [14], valid for high subcooling to 15% quality, is 
recommended. 

Whenever subcooled flow boiling exists to any appreciable extent, 
the cladding surface temperature can no longer be calculated from Newton's 
law of cooling. This is because the heat transfer coefficient in Newton's 
law of cooling does not reflect the increased agitation in the convective 
(now boiling) film due to the growth of the vapor bubbles. Instead, the 
correlation of Jens and Lottes [13] for subcooled boiling, 

- (65) 

may be used. In this relation, AT is the difference between the 

sat 

surface temperature and the local saturation temperature. 
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5.0 The PWR ThermalrHydraullcs Code 

The PWR Thermal-Hydraulics Code calculates the thermal-hydraulics 
performance parameters discussed in Section 4 for the PWR model described 
in Section 3v The .basic operational steps of the code are listed below. 

1. Real«^nteger conversions 

2. Definition of statement functions 

3. ^Accept and print input 

4. Calculate flow area and equivalent diameter of 
coolant channel 

5. Determine calculation increment 

6. Calculate inlet pressure loss 

7. Initialize to 1st calculation increment 

8v .Calculate bulk temperature rise and average 
temperature for 1st increment 

9.- Calculate average coolant properties for 1st 
increment 

10. Calculate pressure loss for 1st increment 

11. Calculate average pressure and saturation 
temperature for 1st increment 

12. { Calculate heat transfer coefficient for 1st 

increment 

■ ^* ■ 

13. Calculate cladding outer surface temperature for 
1st increment 

14. Check for bulk boiling, surface boiling, critical 
heat flux for 1st increment 

15. Calculate temperature distribution in fuel for 
1st increment 

16. Print 1st increment output 

17. Repeat steps 7-16 for remaining increments 

18. Calculate exit pressure loss and total pressure loss 

19. Print coolant exit conditions. 
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The details of these specific steps requiring further explanation are 
found in the discussion below. A flow chart showing the calculations and 
logic is given in Section 5.4. 

Most of the statement functions are used for convenience in 
calculating property data. These statement functions are in the form of 
some property, either of the coolant, cladding, or fuel, as a function 
of the corresponding temperature. All of the function forms are polynomials 
and were obtained by fitting polynomials of various degree to tabulated 
data. The order of these polynomials varies, having been selected to insure 
that it is accurate to within 2% over the range of interest of the para- 
meter. Similar polynomial fit statement functions are used to relate the 
local saturation temperature and the "constants" C and m in Equation 64 to 
the local coolant pressure. 

A detailed discussion of the input requirements is deferred until 
Section 5.2. For the moment it suffices to note that this input is 
supplied on three cards; one each for the geometry, the inlet flow 
conditions, and the reactor power level. 

Since the cross section of the model flow channel does not change 
with axial position, the flow area and equivalent diameter are constant 
axially. The flow area is calculated from 



and the equivalent diameter from 
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where the wetted perimeter is given by 

In the code, the model coolant channel and fuel rod are sliced into 
short segments which are stacked axially to form the proper length core. 
The actual length of these segments is selected as a compromise between 
the very short segments for which the assumption of constant properties 
over the section axially is accurate and the very long segment which 
minimizes the computer time. An increment size equal to 1/100 of the active 
core length is built into the code. To study the effect of increment 
size, or to save computer time, this can readily be changed by altering two 
lines of the code. 

After the increment size has been computed, the coolant conditions 
at the beginning of the first increment are established. The bulk tempera- 
ture and mass velocity at the entrance of the first increment are set 
equal to the core inlet conditions provided in the input. The pressure 
at the entrance of the first calculation increment is obtained by sub- 
tracting the inlet pressure loss, computed from Equation 26 and Eq -"^tion 
34, from the pressure at the inlet of the core. In addition, counters 
required in the code logic are set equal to 1, denoting the first increment. 
The first major calculation in the code is the determination of the 

bulk temperature rise for the first calculation increment. This is calcu- 

II 

lated from Equation 60 with q^ determined from Equation 63. For small 
calculation increments, the average volumetric thermal source strength can 
be accurately approximated by evaluating Equation 1 at -H/2 + AL/2; i.e., 
at the center of the increment. 

ie 
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The average bulk temperature, T„ , for the increment is taken as the 

D 

arithemetic average of the inlet and exit values. All coolant properties 

used in the pressure drop, heat transfer coefficient and critical heat 

flux calculations are assumed to be functions of temperature only and are 

evaluated at the average bulk temperature. 

Both the friction factor and heat transfer coefficient calculations 

require the Reynolds number of the flow which is calculated from Equation 

31. Then f , and Nu , are calculated from Equations 32 and 55, 
clrc clrc ' 

respectively. To these values, the multiplicative correction factors are 

applied to yield f and Nu for the actual reactor flow. Knowing f , the 

pressure drop for the increment is calculated from Equation 26. The heat 

transfer coefficient is backed out of the Nu definition. Equation 52. 

The average outer surface temperature of the cladding for the increment 

is determined from Newton's law of cooling with q", T„ and h available 

^ s B 

from previous calculations for the increment; i.e.i 

- - ^ 

« T + —5— ( 68) 

^8 ^B ^ h • 

At this point, several comparisons are built into the code which 

check for proper reactor operation. The first comparison is of the 

average bulk temperature to the saturation temperature evaluated at the 

average pressure of the increment. T^ < T is required for subcooled 

D sat 

operation. If improper operation is observed an informative message is 

printed and the calculations are terminated. The second comparison is 

T to T . If T > T , then surface boiling may be present. The 
s sac s sat 

possibility of surface boiling for the increment is indicated by a printed 
message. If surface boiling is indicated by the above criterion, then a 
comparison of q^ to the critical heat flux calculated from Equation 64 
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is made. The ratio ^^/^^ = ^HFR (critical heat flux ratio) is computed 
for all segments in which surface boiling is indicated. If CHFR is equal 
to one or greater, the program terminates after printing out q and q . 
If surface boiling exists but burnout does not occur, the cladding surface 
temperature is recalculated from Equation 65. 

After the coolant channel calculations for the increment have been 
performed, the radial temperature distribution in the fuel is calculated 
for the increment. 

First, the temperature at the outer surface of the fuel is found. 

In the fuel rod model, the outer surface of the fuel is at the same 

temperature as the inner surface of -the cladding; therefore, the outer 

surface fuel temperature, T^, can be obtained by adding the temperature 

drop across the cladding to the outer surface temperature of the cladding, 

T , which was calculated earlier, 
s' 

The temperature drop across the cladding is obtained from Equation 

22. In Equation 22 the thermal conductivity of the cladding appears. 

This property is temperature dependent and should be evaluated at the 

mean temperature of the cladding. However the V'^cn cladding temperature 

is dependent on the temperature rise across the f ladding. Thus, an 

iterative-type solution is called for. The procedure used is to 

evaliiate k _ at T , calculate T^, calculate T^- evaluate k 

C-I s f cla stiean c^ii 

at T - ^ , compare k __ to k ... If k ^r within 1.0% of k , 

clad mean' c-II c-I c-^II c-i 

accept value of T-. If k __ is not within 1.0% k , let k = 
*^ f c-II c-I c-I 

k and repeat the sequence. In t-^ls way a adding thermal conductivity 
c— II 

accurate to within 1.0% is used. 

The next step is to calculate the. temperature at the center of the 
fuel from Equation 14. Again a temperature dependent thermal conductivity, 
appears, and an iterative calculation procedure like that described 
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in the previous paragraph is used. The fuel temperature at the quarter, 
half, and three-quarter radii are then calculated from Equation 12. 

After the calculations for the first increment are completed the 
results are printed out. Discussion of the output is deferred until Section 
5.3. It should be noted here, however, that in the output tabulation 
the temperatures are reported as having occurred at the axial center of 
the increment. In other words, for the first increment the results are 
reported at Z = H/2 + L/2. After the results for the first increment 
have been printed out, the code increments to the next segment, and continues 
this process until all 100 segments have been spanned. 

After the completion of the calculations and printout for the 100 
segments of the core, the core exit conditions are calculated and 
printed. A calculation of the exit pressure loss using Equations 26 

and 34 is required. Then P and T are computed and printed, ending 

ex ex 

the computations . 



5.1 Nomenclature 



Analysis Symbol 
A 
A^ 



w 



^1' ^2' ^3' ^4 



P 

CHFR 



D 
D 



e 
f 

^circ 
f/f 



circ 



Code Symbol 



AFL 



CP 

CHFR 

DEQ 

DFU 
DCL 

ERR 

F 

FCIRC 
FRAT 



Description 

Area 

Flow area of coolant channel 

Heat transfer area normal to n- 
direction 

Heat transfer area at radius r 

Heat transfer area at outside 
of cladding 

Surface area of wall 

Constant in Equation 64 

Integration constants 

Constant pressure specific heat 

Critical heat flux ratio 

Diameter 

Equivalent diameter of coolant 
channel 

Outside diameter of nuclear fuel 
Outside diameter of cladding 
Energy release per fission 
Error ratio 
Friction factor 

Friction factor for circular tubes 

Ratio of actual to circular tube 
friction factor 

Dimensional constant = 32.174 
Ibm-ft 

sec^""lbf 
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Analysis Synibol 
G 



Code Symbol 
GIN 



Descrip^tion 
Mass velocity of coolant 



H 
H 

e 

i 
i 
k 



in 



k 



c-I 



c-II 



k. 



K 
L 

e 
AL 

m 

• 

m 
n 
N 



IITC 

H 
HE 



KBC 

KCL 
KCLl 

KCL2 

KFU 
KFUl 

KFU2 



DELL 

M 
MD0T 



ff 



Convective heat transfer 
coefficient 

Height of active core 

Extrapolated height of core 

Enthalpy 

Enthalpy at inlet of core 

Thermal conductivity 

Thermal conductivity of the 
coolant 

Thermal conductivity of cladding 

Thermal conductivity of cladding 
at pre-iteration average cladding 
temperature 

Thermal conductivity of cladding 
at updated average cladding 
temperature 

Thermal conductivity of fuel 

Thermal conductivity of fuel - at 
pre-iteration average fuel 
temperature 

Thermal conductivity of fuel - at 
updated average fuel temperature 

Flow resistance coefficient 

Extrapolation length 

Length increment 

Exponent in Equation 64 

Itass flow rate 

Exponent in Equation 55 

Number of fissionable fuel 
nuclei per unit volume 
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Analysis Symb o 1 

Nu 

Nu ^ 
circ 

Option 

p 
p. 



ex 
Pr 

AP 

AP„ 



Si 
% 
% 



It 

s 

I f 1 



Code Symbol 
NU 
NUCIRC 
0P 

PW 

PAV 

PI 

P2 

PEX 

PIN 

PR 



DELPE 
DELPF 



CHF 
SHFAV 
QTPAV 
QTPO 



Description 
Nussclt number 

Nusselt number for circular tubes 

Option = 1 for cosine. Option = 
2 for coupled program 

Wetted perimeter 

Pressure 

Pressure at beginning of increment 
Pressure at end of increment 
Core exit pressure 
Core inlet pressure 
Prandtl number 
Turbulent Prandtl number 
Pressure drop 

Exit or entrance pressure loss 

Frictional pressure loss 

Heat transfer rate 

Heat transfer rate at 

Heat transfer rate in n-direction 

Heat transfer rate at 

Heat transfer rate at wall 

Critical heat flux 

Heat flux at D 

s 

Volumetric thermal source strength 

Volumetric thermal source strength 
at center of core 

Fuel pin radial coordinate 

Radius of fuel 
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Analysis Symbol 
r 

s 
R 

Re 

S 

S/D 

s 

St 



Code Symbol 



^1/4 



^Bl 



^B2 



ex 



in 



sat 



w 



B 



n 



RE 
S 

SEAT 



T 
TO 

T14 
T12 
T34 
TBAV 
TBI 

TB2 

TF 

TEX 
TIN 
TSAV 
TSATAV 



UB 



y 
z 



53 



Description 

Outer radius of cladding 

Core radial coordinate 

Reynolds number 

Center-to-center spacing 
between fuel rods 

Ratio of fuel rod spacing to 
fuel rod diameter 

Stanton number *^ 

Temperature 

Temperature at center of fuel rod 

Temperature in fuel at r = 1/4 r^ 

Temperature in fuel at r = 1/2 r^ 

Temperature in fuel at r = 3/4 r^ 

Coolant bulk temperature 

Bulk temperature at beginning of 
increment 

Bulk temperature at end of 
increment 

Temperature in fuel at 
Core exit temperature 
Core inlet temperature 
Cladding surface temperature 
Saturation temperature 
Wall temperature 
Velocity 

Average velocity in flow channel 
Coordinate in n-direction 
Coordinate normal to wall 
Axial coordinate 
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Analysis Symbol 
Z 



m 



H 



e 

y 

V 

p 

o 

T 

T 



W 



Code Symbol 
ZC 



VIS 



RH0 



superscript — 



Description 

Axial location of center of 
increment 

Axial location of maximum cladding 
surface temperature 

Thermal diffustvity 

Roughness size 

Eddy diffusivity of heat 

Eddy diffusivity of momentum 

Angular coordinate 

Absolute viscosity of water 

Kinematic viscosity of water 

Density of water 

Average fission microscopic 
cross section 

Shear stress; also time 

Shear stress at wall 

Neutron flux 

Average for calculation increment 
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5.2 Code Input 

The input to the PWR Thermal-Hydraulics Code is intentionally 
extremely simple. Only three cards are required; the first lists the 
geometrical quantities, the second lists the coolant inlet conditions, 
and the third lists the parameters that describe the volumetric thermal 
source strength distribution in the reactor. 

The required input data cards are listed below. The units of each 
parameter and their input format are included. Sample input data cards 
are shown in Figure 2. 

Card 1 — Geometry 

(in.) F 10.4 

D (in.) F 10.4 

s 

S (in.) F 10.4 

H (in.) F 10.4 
Card 2 — Coolant inlet conditions 

G (Ibm/hr-ft^) E 10.3 

T^^ (F) F 10.1 

P^^ (psia) F 10.1 

Card 3 — Power distribution 

q'^' (Btu/hr-ft^) E 10.3 

H (in.) F 10.1 

e 

Option (pure number) I 10 

On Card 3 the axial volumetric thermal source strength distribution 
through the core is specified. For independent operation of the code 
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Figure 2. PWR Thermal- Hydraulics Code Input, 
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t f f 

(Option » 1) a cosine distribution is built into the code and q ^ and 

are the parameters which indicate its level and period. When the PWR 
Thermal-Hydraulics Code is coupled with other codes (Option = 2) the cosine 
distribution is over-iidden by the actual power distribution which is 
supplied by the main program. 
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5.3 Code Output 

Three sets of. information are printed out. Each is printed under 
cooplltely explicit l>eadings such that there should be no interpretation 
difficulties. * 

The first set of output ipformation consists of a listing of the 
input that was supplied to the code. 

The second and main set of output is a tabulation of the temperatures, 
h, AP, and CHFR for each increment in the core. The temperatures printed 
out are f^. T^,^, T^. T^. T^. T^^^. From this data, 

observations of the maximum fuel and cladding temperatures and their 
locations can be made. 

The third set of data reports the core exit pressure and temperature. 
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5.4 Code Flow Chart 



lCj(T), k^(T) 
kgCT), y(T) 
p(T), Cp(T) 
T33,(P), m(P) 
C(P) 



Input: 
Dj, Dg, S, H, 



q'", H , Option 



Functional relations between 
variables supplied to code 
as polynomials* 



Geometry, inlet coolant 
conditions, and power 
distribution supplied 
to code. 



Output: 



Dj, Dg, S, H 

^in^ ^in 

q'", H , Option 
^o • e* 



Printout of geometry, 
inlet coolant conditions) 
and power distribution. 





4 s 




n D 

s 


D » 
e 




i 


AL = 


100 



Calculation of flow channel 
geometrical parameters . 



Selection of calculation increment. 



ERIC 



59 



- 56 - 

0 



Re 



G/P(T^„) 



f ^ = 0.184/Re 
circ 



0.2 



f - f 



circ f 



circ 



ik) 



^^E " ^ 2g^ ' '-^ " 



T ■ 






^Bl 


in 




Z « 

c 


- H/2 + AL/2 


I - 1 












nz^ 




1 r 


cos jj- 






e 


^B2 ^Bl 


+ 


i'" IID^/AmCpCTgP 


^B " (^Bl 


+ 


^B2>/2 




Calculation of permanent pressure 
loss at inlet and resulting pressure 
at start of active core. 



Initializing temperature, 
position, and increment 
counter for 1st increment. 



Calculation of average bulk 
temperature of coolant for 
increment. 
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Re 



G/P (Tg) 



f , -» 0.184/Re' 
circ 



0.2 



"circ f 



circ ^ s 



^2 = ^ - ^^F 



T = T (P) 
sat sat 



Pr = 



Nu , = 0.023 Pr°'^ Re°-^ 
circ 



Nu = Nu 



Nu 



circ Nu 



circ ^ s 



h = 



Nu kg (Tg) 



X = T + =■ 

s B 4 D h 




Output: 
"Net Void 
Fortnation"/ 




Calculation of pressure drop 

and average saturation temperature 

for increment. 



Calculation of heat transfer 
coefficient and average cladding 
surface temperature for increment, 



Check to see if bulk coolant remains 
subcooled. Error message printed 
and calculations terminated if bulk 
coolant reaches saturation. 
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Output : 

"Subcooled Boiling"^ 





C(P) (-^ 

\io 


e) (^sat 


T ^0.22 


■ 


4D 

8 






CHFR - 


s c 








. , « /-fiv0»25 -P/900 
T - T + 1.9 (ci") e 
s sa ^ 




k - k (T ) 

cl c s 



0 



Check to see If subcooled boiling 
is possible. If possibility is 
detected, a message is printed and 
check for subcooled boiling burnout 
is initiated. If no boiling is 
indicated, program proceeds to 
calculation of fuel and cladding 
temperature distributions* 



Critical heat flux check. If 
positive, message and local 
conditions printed out and 
calculations terminated. 



Recalculation of cladding surface 
temperature for subcooled boiling 
on surface. 
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T, = T + „■ ^ In 
Is 8k - 

cl f 



T = (T + Tj/2 
c a f 



k^, = k^ (T ) 
cz c c 




^fl = w 




T = 
o 


T + 




(T, + T^)/2 


k^2 = (^f) 






k - k 




e = 


k,2 





cl . c2 





Calculation of temperature rise 
across cladding. The thermal 
conductivity of the cladding is 
evaluated at the mean cladding 
temperature by an iterative 
technique. 



Calculation of fuel centerline 
temperature. The thermal con- 
ductivity of the fuel is eval- 
uated at a mean fuel temperature 
by an iterative technique. 



EKLC 
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^1/4 






T - 
0 


256 k^j^ 








^1/2 " 


T ■ 

0 


- 64k^3L 






9 i'"D^ 


^3/4 = 


T 

0 


256 kjj^ 




I « I + 1 



Z = Z + AL 
c c 

^Bl B2 

P « P 
1 2 



Evaluation of fuel temperature 
at the quarter, half, and three 
quarter diameter points. 



Printout of thermal-hydraulic 
parameters for increment. 



Check to see if entire core haa 
been calculated 



Initializing for next calculation 
Increment. 
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T = T 
ex B2 



U 



B P(T ) 
ex 



Re = 



U,D p (T ) 
Be ex 

y(T ) 
ex 



f . = 0.18A/Re 
circ 



0.2 



f = f 



circ f 



circ 



p(T ) 

e 



2g. 




Calculation of pressure loss in 
gas plenum and downstream rod 
support region. 



Printout of coolant core exit 
conditions • 
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5.5 Code Examples 

The following two examples demonstrate the use of the PWR Thermal- 
Hydraulics Code for calculating the thermal-hydraulic behavior of a 
typical PWR. 
Example 7. 

A PWR core is composed of 0.422 in. diameter fuel rods 12 ft long 
and spaced 0.563 in. in a square array. The fuel diameter is 0.3664 in. 
The inlet conditions are 543 F and 2250 psia. The average coolant 
velocity in the inlet of the core is 13.02 ft/sec. The average fuel 
rod produces 5.95 kw/ft. Assume an extrapolation length equal to 5% 
of the core height. Determine the magnitudes and locations of the 
maximum fuel temperature and the maximum cladding temperature. 
Solution. 

The mass velocity, extrapolated core height, and peak volumetric 
thermal source strength must first be computed from the data given. 

G - pU - (47.0 Ibm/ft^) (13.02 ft/sec) (3600 sec/hr) 
B 

= 2.2 X 10^ Ibtn/hr-ft"^ 

H =H+2L=H+ 2(0. 05H) = I.IOH 
e e 

= 1.10 (12 ft) = 13.2 ft = 158.4 in. 

q'" = 2P /[H sin(IIH/Z H )] 
^ o ave e f e 

2 

= 2(5.95 kw/ft) (12 ft) (3413 Btu/lcw-hr)/{ (13 .2 ft) (0.3664 in.) 
(ft/12 in.)^ sin[n(12 ft)/2(13.2 ft)]} 
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- 4.0 X lip Btu/hr-ft^. 
From the ou.. " printout of the code (see next five pages) we find 



T , =- 723.3 F at Z - 0.960 ft 

clad max 



T, , = 1677.6 F at Z = 0.180 ft. 

fuel max 
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PWR THERMAL HYDRAULICS CODE 



INPUT DATA 



FUEL O.D. 
i INCH I 
0.3664 



CLAD 0.0. 
( INCH } 
0.4220 



ROD PITCH 
( INCH I 
0.5630 



INLET MASS VELOCITY 
I LB/HR-FT**2 ) 
0.2200E 07 

VOLUMETRIC THERMAL 
( BTU/HR-FT**3 I 
0.4000E 08 



CORE INLET TEMPERATURE 
( F ) 
543.00 

EXTRAPOLATED HEIGHT 
( INCH ) 
158.40 



ACTIVE CORE LENGTH 
I INCH ) 
144.00 

CORE INLET PRESSURE 
( PSIA ) 
2250.00 



OPTION 



END OF INPUT DATA 



AT THE INLET OF THE CORE 



CHANNEL AVE. VELOCITY REYNOLDS # FRICTION FACTOR 
(FT/SEC) 

13.02 413202-1 0.0139 
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ERIC 



NO Z-LOC. 


TO 


Tl/4 


Tl/2 




T3/4 




T CLAD 


T BULK 


T SAT H T COEF 


DELTA 


P CHF R/ 


I -5.940 


663.3 


660.3 


639.0 




606.7 




550.6 


543 .1 


'652.6 5497.1 




4.94 




0.0 


2 -5.820 


686.3 


682.7 


657.4 




621.2 




552.2 


543.3 


652.6 5497.9 




4.94 




0.0 


3 -5.700 


709.6 


705.4 


675.9 




633.8 




553. 7 


543.5 


632.6 5498.8 




4.94 




0.0 


4 -5.580 


733.2 


728.4 


694.7 




646. 5 




555.3 


543.7 


652.6 5499.9 




4.94 




0.0 


5 -5.460 


757.0 


751.6 


713.6 




659.4 




556.9 


544.0 


652.6 5501.0 




4.94 




0.0 


6 -5.340 


781.1 


775.1 


732.8 




672.3 




558. 5 


544.3 


652.6 5502.3 




4.94 




0.0 


7 -5.220 


806.8 


800.1 


753.1 




686.0 




560. 1 


544.6 


652.6 5503.7 




4.94 




0.0 


8 -5.100 


831.7 


824.4 


772.8 




^.99. 3 




561.8 


544.9 


652.6 5505.3 




4.95 




0.0 


9 -4.980 


856.8 


848.8 


792.8 




712.7 




563.4 


545.3 


652.6 5507.0 




4.95 




0.0 


10 -4.860 


882.2 


873.6 


812.8 




726.1 




565. I 


545.7 


652.6 5508.7 




4.95 




0.0 


11 -4.740 


907.8 


898. 5 


833.1 




739.7 




566« 8 


546.2 


652.6 5510.6 




4.95 




0.0 


12 -4.620 


933.5 


923.5 


853.4 




753.3 




568. 5 


546.6 


652.6 5512.7 




4.95 




0.0 


13 -4.500 


959.4 


948.8 


873.9 




767.0 




570. I 


547.1 


652.6 5514.8 




4.96 




0.0 


14 -4.380 


985.5 


974.1 


894.5 




780.7 




571. 8 


547.6 


652.6 5517.0 




4.96 




0.0 


15 -4.260 


1011.6 


999.6 


915.1 




794.4 




573.5 


548.2 


652.6 5519.4 




4.96 




0.0 


16 -4.140 


1037.8 


1025.1 


935.8 




808.2 




57 5. 2 


548.7 


652.6 5521.8 




4.96 




0.0 


17 -4.020 


1064.1 


1050.7 


956.4 




821.9 




577. 0 


549.3 


652.6 5524.4 




4.97 




0.0 


18 -3.900 


1089.5 


1075.4 


976.3 




834.9 




578. 7 


549.9 


652.6 5527.0 




4.97 




0.0 


19 -3.780 


1115.6 


1100.9 


996.9 




848.5 




580.4 


550.5 


652.6 5529.8 




4.97 




0.0 


20 -3.660 


1141.7 


1126.2 


1017.4 




862.0 




582.1 


551.2 


652.6 5532.6 




4.98 




0.0 


21 -3.540 


1167.6 


1151.5 


1037.8 




875.5 




583.7 


551.9 


652.6 5535.5 




4.98 




0.0 


22 -3.420 


1 193.4 


1176.6 


1038.0 




888.8 




585.4 


552.6 


652.6 5538.5 




4.99 




0.0 


23 -3.300 


1223.0 


1205.4 


1081.1 




903.8 




587.1 


553.3 


652.5 5541.6 




4.99 




0.0 






69 











TF 

10 

566.1 
570.5 
574.8 
579.1 
583.4 
587.7 
592.0 
596.3 
600.5 
604.8 
609.0 
613.2 
617.3 
621.5 
625.5 
629.6 
633.6 
636.9 
640.7 
644.5 
648.3 
651.9 
655.5 
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24 -3<,180 


1248.9 


1230-5 1101.4 


917.1 




659.0 


588.8 


554-0 


652-5 5544.8 


4.99 


0.0 




25 -3.060 


1274-5 


1255.4 1121.5 


930.2 




662.5 


590.4 


554-8 


652-5 5548.1 


5.00 


0.0 




26 -2.940 


1299.7 


1280-0 L141.3 


943.2 




665.9 


592.1 


555-5 


652.5 5551.4 


5.00 


0.0 




27 -2.820 


1324-6 


1304.2 1160.8 


955-9 




669.2 


593.7 


556-3 


652.5 5554.8 


5.01 


0.0 




28 -2.700 


1349.1 


1328.0 1179.9 


968.5 




672.4 


595.3 


557.1 


652.5 5558.2 


5.01 


0.0 




29 -2.580 


1373.1 


1351.4 1198.7 


980.7 




675.6 


596.9 


557.9 


652.5 5561.7 


5.01 


0.0 




30 -2o460 


1396.5 


1374-2 1217.1 


992.7 




678.6 


598.5 


558.8 


652.5 5565-3 


5.02 


0.0 




31 -2.340 


1419.4 


1396.4 1234.9 


1004.4 




681.6 


600.0 


559.6 


652-5 5568.9 


S.02 


0.0 




32 -2.220 


144L.6 


1418.0 1252.3 


1015.7 




684. 5 


601.6 


560-5 


652.5 5572.5 


5.03 


0.0 




33 -2.100 


1463.1 


1438-9 1269.1 


1026.7 




687.3 


603.1 


561-4 


652.5 5576-2 


5.03 


0.0 




34 -1.980 


1483.8 


1459-1 1285.4 


1037.3 




690. G 


604.6 


562.3 


652.5 5579-9 


5.04 


0.0 




35 -1.860 


1503.7 


1478.5 1300.9 


1047-5 




692.6 


606. 1 


563.2 


652.5 5533.7 


5.04 


0.0 




36 -1.740 


1522-8 


1497.0 1315.9 


1057.2 




695-2 


607.5 


564-1 


652.5 5587-5 


5.05 


0.0 




37 -1.620 


1540-9 


1514.7 1330.1 


1066.5 




697. £ 


608.9 


565-0 


652.5 5591.3 


5.06 


0.0 




38 -1.500 


1558.1 


1531.4 1343.5 


1075.4 




699.? 


610.3 


565.9 


652.5 5595.1 


5.06 


0.0 




39 -1.380 


1574.3 


1547.1 1356.2 


1083-7 




702.1 


611.6 


566-8 


652-5 5598.9 


5.07 


0.0 




40 -1.260 


1589.4 


1561.8 1368.1 


1091-5 




704,2 


613.0 


567.8 


652.5 5602.7 


5.07 


0.0 




41 -1.140 


1603-4 


1575.5 1379.1 


1098-7 




706.2 


614.3 


568.7 


652-5 5606.6 


5.08 


0.0 




42 -1-020 


1616-3 


1588.0 1389.2 


1105.5 




708.1 


615.5 


569.7 


652.5 5610-4 


5.08 


0.0 




43 -0.900 


1628-0 


1599-5 1398.5 


1111.6 




709. S 


616.8 


570.6 


652-5 5614.3 


5.09 


0.0 




44 -0.780 


1638.6 


1609.7 1406.8 


1117.2 




711. e 


617.9 


571-6 


652.5 5618-1 


5.10 


0.0 




45 -0.660 


1647.9 


1618.8 1414.2 


1122-1 




713.2 


619.1 


572-5 


652.5 5621.9 


5.10 


0.0 




46 -0.540 


1656.0 


1626.7 1420.7 


1126-5 




714.1 


620.2 


573.5 


652.5 5625.7 


5.11 


0.0 




47 -0-420 


1662.8 


1633.4 1426.1 


1130.2 




716-C 


621.3 


574.5 


652-5 5629.5 


5.11 


0.0 




48 -0.300 


1668.4 


1638.8 1430.6 


1133-4 




717-^ 


622.4 


575-4 


652.5 5633-2 


5.12 


0.0 




49 -0-180 


1672.6 


1643.0 1434.1 


1135.9 




718.< 


623.4 


576.4 


652.5 5636.9 


5.13 


0.0 




50 -0.060 


1675.6 


1645.9 1436.6 


1137.7 




719.' 


624.3 


577.3 


652.5 5640-6 


5.13 


0.0 
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51 0.060 


1677. 3 


1647. 5 1 4*^^-0 


1 1 "ao n 




7 OA O 


625.3 


578.3 






U .0 




52 0.180 


1677.6 








7a1 . 1 


626. 1 


579.2 


652. 5 5647. 8 




U . U 


53 0.300 


1676.7 


1647.0 1437-Q 








627.0 


580.2 


652. 5 5651. 4 


5 1 


u. u 


5A 0.^^20 


1674.4 


1644. 8 1436.4 


1 1 A Q 

J. O ^ 7 




722.3 


627.8 


581 .1 


652. 5 5654. 9 




u .u 




55 0.540 


1670.9 


1641.4 1433-8 


1 1 ^7 - 5 




7 99 7 
f AA . r 


62 8. 5 


582.1 


'')52.5 


5 16 


u . u 




56 0.660 


1666.0 


1636. 7 1430.3 


11^5 6 




7 9^^ n 

f A^ . U 


629.3 


583.0 


652. 5 566 1. 8 


5-17 


u . u 




57 0.780 


1659.9 


1630. 8 1425.8 






7 9 d 9 
r A 3 . A 


629.9 


563 .9 


652. 5 ' 365. 1 


5 1 A 


U. U 


5*! 0.900 


1652.6 


iL623.7 14?0-^ 


1 1 PQ Q 




r A3. 3 


630. 5 


584.8 


652.5 5668.4 


? . X O 


u . u 




59 1.020 


1644.0 


1615.3 1413.8 


1 1 9A 1| 
X lilO . u 




r A3. 3 


631. 1 


585.8 


652.5 5671-7 


5 1 Q 


n rt 
U . U 




60 1.140 


1634. 1 


1605- 0 ' ^nfi.4 


1191 7 




r A3 . X 


631. 7 


586.7 


652- 5 5A74- fl 


S PO 


A A 




61 1.260 


1623. 1 


1595. 1 13<?a. 1 


1 1 1 A 7 




7 9 9 f) 
f A A . 0 


632.2 


587.6 


652. 5 5677.9 


5 P O 


u . u 




62 1.380 


1611.0 


1583.4 13afl.Q 


1111 9 
1 1 X X . ^ 




7 9 9 A 
1 AA. «f 


632.6 


588.4 




3 . ^ X 


u .u 




63 1.500 


1597.7 


1570.5 1378.8 


1 ] 05 1 




7 91 Q 
f A X . 7 


633.0 


589.3 


652. 5 5683. 9 


5 PI 


u . u 




64 1.620 


1583.3 


1556-5 1^67-fi 


1 no A A 

X U 7 O 




7 9 1 

r A 1 . 3 


633.4 


590.2 


652- 5 5686. fl 




u. u 




65 1.740 


1567.9 


1541. 5 1356. 1 


1 OQl 




720.6 


633.7 


591 .0 


652. 5 5689. 6 




n A 




66 1.860 


1551.5 


1525.6 1343.5 


^ OA'^ A 




7 1 Q 7 

r 1 7 . r 


633. 9 


591 .8 


652. 5 5692.4 


5 P ^ 


u . u 




67 1.980 


1534. 1 


1508. 7 1 ^*^n 7 


1 07S A 




718.7 


634. 1 


592.7 


652.4 56Q5-.0 


*\ 9 A 


A n 




68 2.100 


1515.7 


1490.9 13x6.2 


1 oaa a 




7 17 iL 

r X r .0 


63 4.3^ 


593.5 


652-4 56Q7-6 


(> Pa 


u. u 


69 2.220 


1496. 5 


1472 .3 1 "^n 1 ^ 


1 OS / 7 




/ Xo. D 


634. 4 


594.3 


652 -4 5700-P 


S 9 S 


A A 




70 2.340 


1476. 5 


1452.8 I>fl6-P 


1 HAA 9 




715.1 


634. 5 


595.0 


652-^ 570^^6 


5 PA 
^ . A O 


u . u 


71 2.460 


14 5.7 


1432. 6 1270.2 


1 O"? ft -'^ 




7 1 'a 7 

r 1 3 . f 


634. 5 


595.8 


652-4 570 5-0 


5-P A 

J . A O 


u . u 




72 2*580 


1434, 2 


1411. 7 1 253. 7 


1 OP A I 

X U A O . JL 




719 9 
r 1 A . A 


634.5 


596.5 


652.4 5707.3 


5-P7 


o o 




73 2.700 


1412.0 


1390. 2 12' 6-6 






7 1 n *\ 


634.4 


597.3 


652.4 5709r 5 


5-27 


0-0 


74 2.820 


1389.2 


1368.1 1219.1 


1 on A - 5 




7nA A 


634.3 


598.0 


652.4 5711.6 


5.28 


0-0 




75 2.940 


1365.9 


1345.4 1201.1 


995.2 




706.9 


634.2 


598.7 


652.4 5713.7 


5.28 


0.0 




76 3.060 


1342.0 


1322.2 1182.8 


983.7 




705.0 


634.0 


599.3 


652.4 5715.6 


5.29 


0.0 




77 3.180 


1317.7 


1298.6 1164.0 


971o9 




702.9 


633.7 


600.0 


652.4 5717.5 


5.29 


0.0 
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78 3.300 


1293.0 


1274.6 1144.9 




959.9 


0.0 


700.8 


633.4 


600.6 


652.4 5719.4 




5.30 




698.5 


79 3.420 


1264.0 


1246.4 1122.7 




945.9 




633. 1 


601.3 


652.4 5721.1 




5.30 




0.0 


696.1 


80 3.540 


1239.2 


1222.3 1103.4 




933.7 


0.0 


632.7 


601.9 


652.4 5722.8 




5.31 




693.7 


81 3.660 


1214.1 


1197.9 1084.0 




921.4 




632.3 


602.4 


652.4 5724.4 




5.31 




0.0 


691.1 


82 3.780 


1188.7 


1173.2 1064.3 




908.8 




631.6 


603.0 


652.4 5725.9 




5.32 




0.0 




83 3.900 


1164.0 


1149.2 1045.2 




896.8 




689.1 


631.3 


603.5 


652.4 5727.4 




5.32 




0.0 




84 4.020 


1138.2 


1124.1 1025.2 




884.0 




686.3 


630.7 


604.0 


652.4 5728.8 




5.33 




0.0 




85 4.140 


1112.3 


1099.0 1005.1 




871.1 




683.4 


630.1 


604.5 


652.4 5730.1 




5.33 




0.0 




86 4. '2 60 


1086.4 


1073.8 


984.9 




858.1 




680.5 


629.5 


605.0 


652.4 5731.3 




5.33 




0.0 




87 4.380 


1060.4 


1048.5 


964.7 




845.0 




677.4 


628.8 


605.5 


652.4 5732.5 




5.34 




0.0 




88 4.500 


1034.5 


1023.3 


944.4 




831.9 




674.3 


628.0 


605.9 


652.4 5733.6 




5.34 




0.0 




89 4.620 


1008.6 


998. 1 


924.2 




818.8 


0.0 


671.1 


627.3 


606.3 


652.4 5734.6 




5.34 






90 4.. 740 


•*982.8 


973.0 


904.0 




805.6 




667.9 


626.5 . 


606.7 


652.4 5735.6 




5»35 




0.0 


664.5 


91 4.860 


» 957^^1 


948.0 


883.9 




792.5 




625.C>>' 


607.0 


652.4 5736.5 




5.35 




0.0 




92 4.'980 


931.5 


923.1 


863.9 




779.4 


0.0 


661.1 


624.8 


607.4 


652.4 5737.3 




5.35 




657.6 


93 5.100 


906.1 


898.3 


844.0 




766.3 


0.0 


623.8 


607.7 


652.4 5738.1 




5.36 




654.1 


94 5.220 


880.8 


873.8 


824.2 




753.3 




622.9 


608.0. 


652.4 5738.8 




5.36 




0.0 


650.5 


95 5.340 


855.8 


849.4 


804.5 




740.3 




621.9 


608.3 


652.4 5739.4 




5.36 




0.0 




96 5.460 


830.1 


824.4 


784.3 




727.0 




646.^ 


620.9 


608.5 


652.4 5740.0 




5.36 




0.0 


643.1 


97 5.580 


805.9 


800.8 


765.2 




714.3 




619. 8 


608.7 


652.4 5740.6 




5.36 




0.0 


639.3 


98 5.700 


781.8 


777.3 


746.2 




701.6 




618.7 


608.9 


652.4 5741.0 




5.37 




0.0 


635.5 


99 5.820 


757.9 


754.1 


727.3 




689.0 


0.0 


617.6 


609.1 


652.4 5741.4 




5.37 






100 5«940 


734.3 


731.1 


708.6 




676.5 


0.0 


631. e 


616.5 


609.2 


652.4 5741.8 




5.37 






EXIT 


PRESSURE , 


EXIT TEMPERATURE 











2244.8 609.3 



POWER PER UNIT LENGTH CKW/FT) 5.95 

O 
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Example 8, 

The effect of oversized cladding on PWR thermal-hydraulic performance 
is to be studied. Consider the PWR described in Example 7 but with an 
outer cladding diameter of 0.500 in. (instead of 0.422 in,). Compare the 
following thermal -hydraulic parameters for the "oversized cladding" and 
"normal" flow channels: 

a. Location and magnitude of maximum fuel temperature 

b. Location and magnitude of maximum cladding temperature 

c. Exit water temperature. 
Solution. 

Each flow channel in the fuel assembly wii: pxpnrience the same 
pressure drop. The "oversized cladding" flow channel will have a smaller 
equivalent diameter and hence will pass less flow for the same pressure 
loss (as shown in Example 5) . To assure approximately the same pressure 
loss as in the "normal" channels the mass velocity of the "oversized 
cladding" channel must be adjusted. The equivalent diameters are computed 
first. 

"Normal" channel: 

= (0.563)^ - |- (0.422)^ 
= 0.1771 in.^ 



p = no = n(0.422) 
•^w s 



= 1.327 in. 



D = 4 A./p 
e f w 



79 



- 4 (0.1771)/1.327 

- 0.534 in. 



"Oversized cladding" channel: 



«2 n ^2 



(0.563)^ - f (0.500)2 



0.1206 in.^ 



p - HD - n(0.500) 

s 



= 1.571 in. 



D - 4 A./D 
e f e 



-= 4(0.1206)/1.571 
« 0.307 in. 
From the results of Example 5, 

v2/3 

G » G„(D /D „) 
oc N eoc eN 

= 2.2 X 10^ (0.307/0.534)2^^ 



6 2 
= 1.52 X 10 Ibm/hr-it . 



The input to the PWR Thermal-Hydraulics Code is the same as that 
for Example 7 except for the new value of G computed above and the 
oversized D . From the code (printout next seven pages) the following 



s 

results are obtained: 
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PWR THEftHAL HYDRAULICS CODE 
INPUT DATA 



FUEL 0.0. 
< INCH ) 
0.3664 



CLAD U.D. 
( INCH ) 
0.5000 



ROD PITCH 
( INCH I , 
0.5630 



INLET MASS VELOCITY 
I LB/HR-FT**2 ) 
0.1520E 07 

VOLUMETRIC THERMAL 
I BTU/HR-FT**3 ) 
0.4000E 08 



CORE INLET TEMPERATURE 
( F ) 
543.00 

EXTRAPOLATED HEIGHT 
( INCH ) 
158.40 



ACTIVE CORE LENGTH 
( INCH ) 
144.00 

CORE INLET PRESSURE 
f PSIA ) 
2250.00 



OPTION 



END OF INPUT DATA 

« 4i :^ 4i « 4i « 4i 4i4t «c 4i« 4i « 4i 4i 4i « 4i 4i « 



AT THE INLET OF THE CORE 



CHANNEL AVE. VELOCITY REYNOLDS # FRICTION FACTOR 
CFT/SEC) 

8.99 164104.4 0.0167 
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NO Z-LOC. TO 
T CLAD T BULK 



Tl/4 Tl/2 
T SAT H T COEf 



T3/4 TF 
DELTA P CHF RATIO 



1 -5.940 


684.9 


661. 8 


660.2 








551.9 


543.2 


652.6 3993«5 




4»31 




U « U 


2 -5.820 


712.3 


708.7 


682.8 




64d. 8 




553.9 


543.6 


652.6 3994.7 




4.31 




O.U 


3 -5.700 


740 .2 


735.9 


705.7 




662.6 




555.9 


544.0 


652.6 3996. 1 




4.31 




0.0 


4 -5.580 


768.5 


763.6 


729.0 




679 .!> 




558.0 


544.5 


652.6 3997.7 




4.32 




0.0 


5 -5.460 


796.6 


791.0 


751.9 




696«U 




560. L 


545 .1 


652.6 3999*5 




4.3Z 




U .0 


6 -5.340 


825. 5 


dl9.3 






7io . 1 




562.2 


545 .7 


652.6 4001. 5 




4.3Z 




U .U 


7 -5.220 


856.4 


849.4 


800.7 




731 .1 




564.4 


546.4 


652.6 4003. 7 




4.32 




u.o 


8 -5.100 


886.4 


878.8 


825.2 




746.6 




566.6 


547.1 


652.6 4006. 1 




4. 33 




0«0 


9 -4.980 


916.9 


908.6 


850.0 




766.4 




568.9 


547.9 


652.6 4008.6 




4.33 




0.0 


10 -4.360 


947.7 


938.7 


875.1 




784.2 




571.2 


548^8 


652.6 40X1.3 




4.34 




0.0 


11 -4.740 


978.9 


969.1 


900.4 




802.2 




573.6 


549.7 


652o6 4014.2 




4.34 




0.0 


12 -4.620 


1010.4 


999.9 


925.9 




820. 3 




576.0 


550.7 


b52«6 4017.3 




4.34 




CO 


13 -4,500 


1042. 1 


1030.9 


951.6 




838.5 




578.4 


551.7 


652.6 4020*5 




4.35 




0.0 


14 -4.380 


1074.2 


1062*i 


977.5 




856.7 




580.8 


552.8 


6^)2.6 4023.9 




4.35 


675.0 


0.0 


15 -4 ,d60 


1106.4 


1093.6 


?.003.6 






583.? 


553.9 


652.6 4027.4 




4.36 




O.D 


16 - t. 140 


1138.8 


1125.2 


1029.7 




893.3 




58"^. 7 


555. L 


652.6 4031.1 




4.36 






17 -4.020 


1171.4 


1157.0 


1055.9 




911.7 




£j8.2 


556.3 


652.6 4034.9 




4o37 




0.0 


r -3.900 


1204.0 


1188.8 


1082.2 




930.0 




590.8 


557.6 


652.6 4038.9 




4.38 




0.0 


i9 -3.780 


1236.6 


1220.7 


iioe.5 




948.3 




593.3 


558.9 


652.6 4042.9 




4.38 




0.0 


.0 -3.660 


1269.3 


1252.5 


1134.7 




966.5 




595o8 


560.2 


652.6 4047.1 




4.39 




0.0 


r:- -3.540 


1301.8 


1284.3 


1160.8 




984.6 




>9<3»4 


> ' o 


652o6 4051«4 




4.40 




OeO 


-3.42.C 


1334.3 


1315.9 


1186.8 


1002.5 




600.9 


563.1 


652.6 4055.1 




4.40 




0.0 



586.1 , 

594. 1 ' 

602.2 I 

610.3 ' 
617*b] 
625. T 
633.6 

6'Vl.J? j 

I 

649r- ?i 
?.5 >'.i ' 

664. a 

6 » it • 5 
680.1 
687.6 
695 o \ 
752.4 
709.7 
716.9 
724.0 
730-9 
737.8 
744.5 



7G 
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X*f HO 93 L£.^D mD 


XU r P • P 
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U. U 
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1 47A Q 1 '31 1 A A 


1 no o A 
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r r o. 3 
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65P A 4n7A A 


A A A 


n n 
u .u 




Pft — ?-.7on 


1 53n p 


1 5n7 O 1 '^14^ *? 
LD\J 9 mU LD*rDmC, 


1 1 no K 
X XU 7 • P 




7 A 9 O 

r o A«4C 


616.3 


57P 4 

D § £. m^ 


A5P a A.nA'^ A 
03^*0 HU03«0 


A A <; 


n A 
U .U 




7 ^ m J 0\J 


1 5A0 5 


1 (>'^A 4 1 3A7 "a 






roo. U 


61 R. fi 




A5P A 4nfiA ^ 


A A A 


n n 
U .U 






1 5QO O 


1 *>Afi P 1 '^on C 
X30P«^ X37U«V 


1 1 A o n 

1 IH^ m\J 




7Q^ a 
#73 .O 


6P 1-4 


575 - A 


A5P-A 4nQ^ 1 


A A 7 


n n 
U . U 




31 — P-34n 


161 A- 7 

■X O X O • f 


1 50*^ P 1 ^ A 
X:973«^ 1H13«0 


1 1 *i7 7 

X X D r • 1 




r 77.U 
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A5P A 40Q7 Q 
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4 4 A 
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3P — P- PPO 


1 A4A 6 


1 APn 4 1 4 3 A n 


1 1 79 A 
XX 9 CmO 
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OUH. D 
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57Q 3 

■ 7 O 3 


A5P A 41 HP 7 


A AO 
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u .u 




33 ~ p 1 nn 


1 A7 3 5 
XO 9 D mD 


1 A4A 7 1 Ac;7 c; 
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4.58 
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42 -1.020 
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649.6 


597.4 
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4.61 


0.0 




43 -0.900 


1678.7 


1846.6 1621.8 


1300.6 




850.9 


651.9 


599.3 


652.5 4151.8 


4.62 
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SUBCOOLED BOILING 








44 -0.780 


1694.8 


1862.5 1635.0 


1310.3 




655.8 


655.9 


601.1 


652.5 4015.2 


4.63 


0.18 





SUBCOOLED BOILING 
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45 -0.O60 1903.2 1870.6 1641.6 1314.5 856.7 
655.9 603.0 652.5 4176.3 4.64 0.18 

SUBCOOLEO BOILING 

46 -0.540 1910.2 1877.5 1647.1 1318.1 857.5 
655.9 604.8 652.5 4345.4 4.66 0.18 

SUBCUOLEO BOILING 

47 -0.420 1915.9 1883.0 1651.4 1320.9 858.1 
655.9 606.7 652.3 4523.1 4.67 0.18 

SUBCOOLEO BOILING 

48 -0.300 1920«l 1887.1 1654.7 1323.0 858.6 
655.9 608.5 652.5 4710.3 4.68 0.19 

SUBCOOLEO BOILING ! 

I 

49 -0.180 1922.9 1889.8 1656.9 1324.4 858.9 
655.9 610.4 652.5 4907.8 4.70 0.19 

SUBCOOLEO BOILING 

50 -0.060 1924.3 1891.2 1658.0 1325.1 859.1 
655.9 612.2 652.5 5116.7 4.71 0.19 

SUBCOOLEO BOILING 

51 0.060 1924.3 1891.2 1658.0 1325.1 859.1 
655.9 614.0 652.5 5338.1 4.72 0.19 

SUBCOOLEO BOILING 

52 0.180 1922.9 1889.8 1656.9 1324.4 858.9 
655.9 615.8 652.5 5573.5 4.74 0.19 

SUBCOOLEO BOILING 

53 0.300 1920.1 1887.1 1654.7 1323.0 858.6 
655.9 617.6 652.5 5824.1 4.75 0.20 

SUBCOOLEO BOILING 

54 0.420 1915.8 1882.9 1651.4 1320.9 858.1 
655.9 619.4 652.5 6092.1 4.76 0.20 

SUBCOOLEO BOILING 
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55 0.540 1910.2 1877.5 1647.0 1318.1 857.5 
655.9 621.1 652.5 6379.3 4.78 0.20 

SU8C00LE0 BOILING 

56 0.660 1903.2 1870.6 1641.5 1314.5 856.7 
655.9 622.8 o52.5 6688.2 4.79 0.20 

SUBCOOLEO BOILING 

57 0.780 1894.7 1862.4 1635.0 1310.3 855.7 
655.9 624.6 652.5 7021..8 4.81 0.20 

SUBCOOLEO BOILING 

58 0.900 1884.9 1852.9 1627.4 1305.4 854.6 
655.9 626.3 652.5 7383.5 4.t52 0.21 

SUBCOOLEO BOILING 

59 1.020 1873.8 1842.0 1618.7 1299.8 853.3 
655.9 627.9 652.5 7777.3 4.83 0.21 

SUBCOOLEO BOILING 

60 1.140 1861.3 1829.9 1608.9 1293.5 851.9 
655.9 629.6 652.5 8208.2 4.85 0.21 

SUBCOOLEO BOILING 

61 1.260 1847.5 1816.4 1598.2 1286.6 850.3 
655.9 631.2 652.5 8682.7 4.86 0.21 

SUBCOOLEO BOILING 

62 1.380 1832.3 1801.7 1586.4 1279.0 848.6 
655.8 632.8 652.5 9207.8 4.8? 0.21 

SUBCOOLEO BOILING 

63 1^500 1816.0 1785.8 1573.6 1270.7 846.6 
655.8 634,4 652.5 9793.3 4.89 0.21 

SUBCOOLEO BOILING 

64 1.620 1798.3 1768.7 1559.9 1261.9 844.6 
655.8 636.0 652.5 10451.2 4.90 0.22 

SUBCOOLEO BOILING 
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65 1.740 1779.5 1750.4 1545.2 1252.4 842.4 
655.8 637.5 652.5 11196.7 4.92 0.22 | 

SUBCOOLEO BOILING 

66 1.860 1759.5 1730.9 1529.7 1242.3 840. ol 
655.8 639.0 652.5 12050.6 4.93 0.22 ■ 

SUBCOOLEO BOILING | 

67 1.980 1738.4 1710.4 1513.2 1231.7 837.5 
655.8 640.5 652.5 13039.4 4.94 0.22 p 

SUBCOOLEO BOILING 

68 2.100 1716.3 1688.9 1495.9 1220.5 834.91 
655.8 641.9 652.5 14200.1 4.96 0.23 ■ 

SUBCOOLEO BOILING | 

69 2.220 1693.1 1666.3 1477.8 1208.8 832.1 
655.8 643.4 652.5 15585.5 4.97 0.23 ^ 

SUBCOOLEO BOILING 

70 2.340 1668.9 1642.8 1459.0 1196.6 829. iH 
655.7 644.7 652.5 17270.8 4.98 0.23 ■ 

SUBCOOLEO BOILING | 

71 2.460 1643.9 1618.4 1439.4 1183.9 826.1 
655.7 646.1 652.5 19370.3 5.00 0.24 _ 

SUBCOOLEO BOILING * 

72 2.580 1618.0 1593.2 1419.2 1170.7 822. 8« 
655.7 647.4 652.5 22064.8 5.01 0.25 | 

SUBCOOLEO BOILING ■ 

73 2.700 1591.3 1567.3 1398.3 1157.1 819.5* 
655.7 648.7 652.5 25658.2 5.02 0.26 ^ 
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74 2.820 1563.9 1540.6 1376.9 1143.2 816.0 
655.7 650.0 652.5 30707.1 5.03 0.28 

SUBCOOLEO BOILING 

75 2.940 1535.8 1513.3 1355.0 1128.9 812.4 
655.6 651.2 652.5 38341.8 5.05 0.31 

NET VOID FORMATION 
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Subcooled (surface) boiling begins in the coolant channel 
at Z - -0.780 ft; i.e., about 3/4 foo*: upstream of the center 
of the core. No boiling was observed for the normal channel. 
The coolant bulk temperature exceeds the local saturation 
temperature at about 3 ft past the center of the core. The 
calculations were terminated at this point since this is an 
improper operating condition for the PWR. The reason for the 
observed behavior is that the "oversized cladding" coolant 
channel is flow starved and can not remove the energy supplied 
to it by the fuel without boiling. 

The maximum fuel temperature is 1924.3 F and occurs at the 
center of the core. This compares to the 1677.6 F at Z = 
0.180 ft for normal core operation. 

The maximum cladding temperature observed was 859.1 F and also 
occurred near the center of the core. This compares to 723.3 
at Z = 0,960 ft for normal operation. 
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5.6 Listing of Code 



C 
C 
C 

C VPI PWR THERMAL HYDRAULICS CGDE. 

C 
C 

ir^TEGE» OP 

REAL KFl!,KCL,<BC,M,rJUCIRC,NU,f<DCT,KCLl,KCL2,KFUl,KF'J? 

C 

C TEKPEPATURE DEPENDENT PROPFRTIPS. 

C 

KFIJ(T)=C.61C21E1-0.^6365E-2«T+0.13063E-5*T**2 
KCL(T)=C.536^3E1+0.2321^E-2*T 

K8C(T )=C. 11711+0. 13910E-2«T-0. 18 102E-5«T**2 
VIS ( T )=0. 85067-0. 17501E-2*T+0.11420E-5*T**2 
PHOCT )=0.5778E2+0. 28018 E-1*T-0.88346E-A*T**2 
CP{ T)=0.47088E1-0.15753E-1«T+0.17233E-4«T*«2 

C 

C SATUPATION TEMPERATURE AS FUNCTION OF PRESSURE. 

C 

TSAT( P)=0.45602E3+0.1l033*P-0.10l99F-4*o**2 

C 

C PRESSURE DEPENDENT •'CONSTANTS'' FCR EQUATION N-64. 

C 

M( P»=-0.^»8 141E- 1+0. 22471 E-3*P+0.83018E-9*P**2 

C(PJ=+0.113l2Ei-0.78004E-3*P+0.33119E-6*P**2-0.56355E-lO*P**3 

c 

C READING INPUT DATA. 

C 

READ(5,5)0FU,DCL,S,H 
5 FQRMAT(3F10.4,F10.1 ) 
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REaD( 'i, 10)GIN,TIN,PIN 
10 F0RMAT(E10.3,2F10.1) 

REA0(5,15)QTP0,HE,CP . 
15 F0RMAT(E10.3,F10.1, 110) 

P'JTNTOUT OF INPUT DATA. 

WRITE (6, 20 J 

20 FORMAT ( IHl f iXf * *'^***********«***************» , /f 3X, » PW^ THcRMAL HY 
<DP AUL ICS CCDE» , / ,2Xf » ****«***i:<«**«******<t ********* f , //,3X, • I NPUT D 

<ATA',//) 
wRITF (6,25)0FUf DCL,S,H 
25 F0R*U.T(3X, 'FUEL C .D . » , 6X , » CLAD O.D, • , 6X, » PCD P I TCH» , 6X , • ACT I VE COR 
<F • MHS/fSXfM INCH )»,7X,M I MCH )«,7X,M INCH )«,12X,M INCH 
< /^4X,F6„4,9X,F6.4,9X,F6.4,L4X,F6.2,/) 
V ^» ''vSOJGIN.TlNf PIN 
3J F .. ..rSX^MNLET MAS^ VFLQC I TY» , 3X , » CORE INLET TEMPERATURE », 3X, • CO 
<RE INLET PRESSURE' ,/f5X,» ( LB/HR-FT**2 )»,13X,M F )«,15X,M PSI 
<A ) »i/f7X,E11.4fl4X,F7.2,17X,r7.2,/) 
WRITE(6,35)QTP0f HE 
35 F0RMAT(3X, »VQLUHETRIC THERMAL' ,4X, 'EXTRAPOLATED HE IGHT ' , / , 4X, ' ( PiT 
<U/HR-FT**3 )',11X,»( INCH ) » , / , 6X , E 1 1 . 4 , 1 5X, F6 ,2 , // ) 
WRITEt6f40)0P 

40 rORMAT(19Xf 'OPTION' ,3X, 12,///, 3Xf 'E NO OF INPUT DATA', 
</ , 3X, ' *****♦♦*♦ ******************♦«*! ,////) , 



CONVERTING INCHES TO FEET. 

DFU=DFU/12.0 
S=S/12.0 

DCL = 0CL/12.0 Q4 
H=H/I2.0 
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HE=HE/12.0 
PI=3. 14159 

C 

C GEOMETRY CALCULATIONS. 

C 

AFL=S**2.0-PI*DCL**2.0/4.0 

PW=PI*OCL 

DEQ=4.0*AFL/PW 

SRAT=S/DCL 

C 

C RATIO CF ACTUAL FRICTION FACTOR OF COOLANT CHANNEL TO CIRCU 

C TUBE FRICTION FACTOR. 

C 

FRAT=-3. 4751+8. 0528*SRAT-4.7047*SRAT**2+0.91621*SRAT**3 

C 

C MASS FLOW PATE CALCULATION. 

C 

MOOT=GIN*AFL 

C 

C SELECTION DF AXIAL CALCULATION INCREMENT. 

C 

0ELL=H/100.0 

C 

C CALCULATION OF PRESSURE LOSS AT INLET OF CORE. 

C 



UB=GlN/RHC{TIN)/3600.0 

RE=UB*DEQ*RHO(TIN)/VIS(TIN)*360C.O 

FCIRC=0.184/RE**0.2 

F=FCIRC*FRAT 

OELPF=F*DELL*RHO(TIN)«UB**2/DEQ/2.0/32.2 
0ELPE=0.5*RH0( TIN )*UB**2/2. 0/32.2 
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CALCULATION OF CONDITIONS AT BEGINNIMG OF FIRST IMCREME'^T. 

Pl=PIN-(OELPF-»-DELPE} 

TB1=TIN 

! = 1 

CriFR=G.O 

CALCULATION OF AXIAL COORDINATE OF CENTER OF FIRST INCREMENT. 
ZC=-H/2.0+DELL/2.0 

CALCULATION OF AVERAGE BULK TEMPERATURE FOR INCREMENT. 
WRITE(6,42) 

42 F0RV,AT(5X,« AT THE INLET OF THE r.ORE»f//> 

WRITE(6,45)U8,RE,FCIRC 
45 FORMATdX, 'CHANNEL AVE. VELCCiTY 2X, 'REYNOLDS # ' , 3X , ' FRI CT 1 ON FA 
<CTOR' f/f7X,» (FT/SEC)' f /,8Xf F6.2,12XfF8.1f8XfF7.4f//» 

WRITE(6«50) 

50 FORMAT ( ' *********«*«*«*4(**««***«*«4i*«****««**4>***4^«4c*****«4^«««**«* 

<♦***• f///) 
KRITE(6f55) 

55 FORMAT <3X, 'NO' f 3X , ' 2-LOC. » i 6X, 'T0',7Xf 'T 1/4' f 6X, ' T 1/2' i 6X, ' T3/4' t 7 
<X,'TF' ,6X, 'T CLAD',4Xf'T BULK'fSX.'T SAT',3Xf»H T COEF 4X ,' DELTA 
<P» ,3X, 'CHF RATIO' ,/) 
60 QTPAV=QTPG*COS(PI*ZC/HE) 

TB2=TBl+QTPAV*Pi*DFU**2*DELL/4.0/MD0T/CP(TBl) 
T8AV=(TBl+TB2)/2.0 

CALCULATION OF PRESSURE DROP AND AVERAGE PRESSURE FOR INCREMENT 
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UB=GIN/PHC(TBAV)/3600.0 
RE=UB*nEQ*RHG(TBAV) /VIS (TBAV) ♦3600,0 
FCIRC=0,184/RF**0.2 
F=FCIRC*FRAT 

DELPF = F*DELL*RH0(TBAV)*UB**2.0/DC-Q/2.0/32.2 

P2=P1-DELPF/144,0 

PAV=(Pl+P2)/2.0 

^ CALCULATIOM OF AVERAGE SATURATION TEHPEkATURE FOR INCREMENT. 

TSATAV=TSAT(PAV) 

CALCULATION OF HEAT TRANSFER COEFFICIENT FOR INCREMENT, 

PR=CP(TBAV)*VIS(TBAV)/KBC{TBAVI 
NUCIRC=0,023*PR**0,4*RE**0,8 
NU=NUCIRC*FRAT 
HTC=NU*KBC(TBAV)/DEQ 

CALCULATION OF TEMPERATURE AT OUTSIDE SURFACE OF CLADDING FOR 
INCREMENT. 

TSAV=TBAV+QTPAV*DFU**2.0/4.0/DCL/HTr. 

CHECK FOR POSSIBILITY OF NET VOID FORMATION. IF DETECTED, "Nt^ 
VOID FORMATION" PRINTED GUT AND CALCULATIONS TERMINATED. 

IF( TB2-TSATAV) 75, 65,65 
65 WRITE(6,70| 

70 FORMAT(/,I0X, • NET VOID FORMAT ! CN« ,// » 
GO TO 180 
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CHECK FDR POSSIBILITY OF SURFACF 3DILING. IF DETECTED, "SUB- 
CnCLFD BOILING" PRINTED OUT AND BOILING BURNOJT CHECK niTIATEO 

75 IF( TSAV-TSATAV)115,80,80 
30 WRITE(S,35) 

85. F0RMAT(/,10X, • SUPCOOL ED BLILING',/) 

CALCULATION CF CRITICAL HEAT FLUX RATIO. 

CHF=C(PAV)*(GTN/1C.0**6)**M( PAV)*( TSATAV-TBAV)*«0.22*10, 0**6 

SHFAV = 0TPAV*0FU«*?.ij/4.':/0CL 

CHFR=SHFAV/CHF 

CHFCK FOR HOI LING BL'RNOUT. IF DETECTED, '•SUBCOOLED BU^MOUT" 
pi^IMTEI) OUT AND CALCULATIONS TERMINATED. 

IF( SHFAV-CHF)110,90,90 
90 WRITE(6,95) 

95 F0RMAT(/,10X, • SUBCOOLED BURNOUT*,/) 

WRITE(6,100) 
100 FORMAT {10X,«CHF« ,20X,«SHF« ,/) 

••RITE(6,L05)CHF, SHFaV 
10 5 rORMAT(6X,E10.3,6X,E10.3 » 

GO TO IfiO;^' 

RECALCULATION OF TEMPERATURE AT OUTSIDE OF SURFACE CLADDING AhiO 
HEAT TRANSFER COEFFICIENT FOR SUBCOCLEO BOILING CASES. 

IIO TSAV=TSATAV+1.9*SHFAV**0.25*EXP(-PAV/900.0) 
HTC=SHFAV/(TSAV-TBAV) 

CALCULATION OF TEMPERATUPF AT CUTSIDE SURFACE OF FUEL FOR 1^- - 
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CREMFMT, DETERMINED BY AOD1..G TEMPERATURE RISE ACROSS CLADDIMG 
TO TPMPEFATURF AT OUTSIDE SURFACE OF CLADDING. AN ITERATION 
TECHNIQUE IS USED TO OBTAiN THE CLADDING THERMAL CONOUr.TIVITV 
TO WITHIN IS. 

115 KCL1=KCL(TSAV) 

120 TF=TSAV+QTPAV*DFU**2.0/8.0/KCLl*Al.OGlOCL/nFU) 

TCLAV=(TSAV+TFl/2.0 

KCL2=KCL(TCLAV> 

£RR=A6Sl (KCL1-KCL2)/KCL2) 

IF(eRP-0«.010)130,130,125 
P.25 KCL1 = KCL2 

GO TO 120 

CALCULATIUM OF TEMPERATURE AT CENTER CiF FUEL FOR INCREMENT. AN 
ITERATIOri TECHNIQUE IS USED TC OBTAIN FUEL THERMAL C3NDUCTIVITV 
WITHIN 1%, 

L^C KF;ji=kFU{TF) 

135 Tu=TF+QTPAV«DFU**2.0/16.0/KFUl 
TFAV=(TF+TO)/2«0 
KFU2=KFU(TFAV) 
FRR=ASS( (KFU1-KFU2J/2.C) 
IF(ERR^0.010)I45*145,1^0 

i^O kf;ji = kfu2 

GO TG 135 

CALCULATION OF FUEL TEMPERATURE AT 1/4, 1/2, AND 3/4 FU^L RADII 
FOP INCREMENT- 

145 T14=Tb-QTPAy*DFU*')'2. 2/16. 0/KFUl/ 16.0 
T12=iT0-yTPAV*PFU**2.0/16.0/KFUl/4.0 
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T3^ = Tn-^/TPA^y'<';FU*«2.0/16. O/K rUl/lo.UvQ.O 

peiNTGUT OF RESULT FCK CALCULATICN INCPEHENT. 

'rtf'.IT-lG, 150)1 , ZC.»TPtTl't,T12vT34,TF,Tf,AV,TBAV,TSATAV,HTC,Dr-LPF,CHFR 
150 FORMAT ( IX, ri,FlC.3»9F10.1f/^X,F6,2,6X,F6.2» 

CHtr.K TO SFF IF CALCULATIDNS HAVE BEEN COMPLETED FOR EMtIRE 
CnPE. 

IF( I-lOO) 15'3» i6c-, 160 

SETTING Co-DlTICr.S FTr: 6EGINMNG Uf NEXT CALCULATKM IVCREME'^T. 

155 1=1+1 

ZC = 7C+'f LL 
icll = T32 
P1 = P2 

TRAMSFtP?!\G TC START CALCULAT lOr^S FOP NEXT INCREMENT. 

CU T". 6C 

C&lCULATI^.^ f f EXIT PRESSURE LOSS . 

1 hO \.JP = 'J: I ri/RHT ( T>:i2 ) /3600.0 

RF. = Li3»-0EQ*KSO(TH?)/VIS(T32)*360C.C 

FCIRC=0.1f* V/Rf **0.2 

F=FCTRC*FRAV 

DELPF=F*0ELL*RH0(TB2)*UB**2/DE0/2. 0/32. 2*6. C 
rA-LP = -t<H0«TP2)*'je**2/2.0/32.2 
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CALrULATICM OF CORE EXIT rGNDITIONS OF COOLANT. 

PEX=P2-( OELPF+DELPE)/ 144.0 
TEX=TB2 

PRINTOUT OF CORE EXIT COOLANT CONDITIONS, 
WRITE (6, 165) 

65 Fl3RMAT(//,10X,«EXIT PR E SSi IP E • , ICX , • EX I T TEMPERATURE*, 

wRITF(6,I70)PEX,TEX 
70 F0RMAT(13X,F6.1,18X,F6.H 

CALCULATION OF AVERAGE POWER OF FUEL ROD IN KH/FT. 
PRINTOUT OF AVERAGE POWER. 

P0VL=0TP0*HE*DFU*»2/2.C*SIN{PI*H/2.G/HE)/H 
P0VL=P0VL/3412 
rtRITE(6, 175)PDVL 
75 F0RMAT(//,5X,« POWER PER UNIT LENGTH (KW/FTI '.Fd.Z, 

60 STOP 

END 
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7.0 Problems 



1. By performing an energy balance on a differential element with sides 
dx, dy, aiid dz, derive the steady state heat conduction equation In 
Cartesian coordinates. 

2. Derive the steady state heat conduction equation in cylindrical 
coordinates. 

3. Develop an expression for the temperature distribution in a long, thin, 
hollow cylindrical fuel element with inner radius r^, outer radius r^, 
and volumetric thermal source strength q"'» Assume that the inner 
surface is Insulated and at temperature T^. 

4. A long, thin, hollow cylindrical fuel element with voltimetric thermal 
sotirce strength q'" is cooled on both the inner and outer surfaces 

r^^ and r = r surfaces) such that the surface temperatures are 
Tj and T^, respectively. Derive an expression for the location of 
the maximum temperature in the fuel. 

5. Calculate the ratio of peak power to avarage power for a fuel rod of 
constant cross section assuming that the volumetric thermal source 
strength varies as q*" « ^o'*' (IIZ/Hg) and that the extrapolation 
length is 10% of the core height. Compare this ratio to that for the 
case where the extrapolation length is assumed to be zero* 

6. Show that for a cylindrical shell without Internal heat generation the 
use of an average constant thermal conductivity evaluated at the 
arithmetic mean of the Ivmer and outer surface temperatures is 
equivalent to a linear variation of thermal conductivity with 
temperature. 

7. A FWR fuel rod consists of 0.500 in. diameter fuel pellets in a 
0.500 in. ID, 0.600 in. OD tube of Zlrcaloy 2 cladding. The heat 
transfer coefficient at the outer surface of the fuel rod is 5000 
Btu/hr~ft -F. The thermal conductivities of the fuel and cladding 
are 1.95 Btu/hr-ft-F and 8.0 Btu/hr-ft-F, respectively. The coolant 
bulk temperature is 550 F and the volumetric thermal source strength 
is 1.2 X 10^ Btu/hr-ft^. Determine the maximum fuel temperature and 
the maximum cladding temperature. 

8. Same as Problem 7 except the ID of the cladding is 0.540 in. The 
0.040 in. diametrical gap between the fuel and cladding is filled with 
helium gas. Assuming that the heat transfer through the helium is by 
molecular conduction only (k^^ = 0.140 Btu/hr-^f t-F) determine the 
maximum fuel and cladding temperature. 

9. A ten foot long PWR fuel assembly contains 0.520 in. OD fuel rods 
spaced 0.680 in. in a square array. The average conditions in the 
.core areg2500 psla and 550 F. The mass velocity of the water is 
2.5 X 10 lbm/hr~f t^. Determine the pressure drop through the- core. 
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10. Determine the effect on pressure drop of a 10% increase in cladding 
OD for the fuel in Problem 9. Assume the spacing, coolant conditions 
and mass velocity remain the same. 

11. Determine the heat transfer coefficient for the fuel and flow conditions 
of Problem 9. 

12. A PWR fuel consists of 0.50 in. diameter UO^ fuel pellets contained 
in 0.50 in. ID, 0.57 in. OD Zirconium 2 cladding. The fuel rods are 
spaced 0.70 in. in a square array. At a particular location in the 
core the pressure is 2500 psia, the temperature is 520 F and the 
coolant velocity is 16 ft/sec. The volumetric thermal source strength 
is 6 X 10^ Btu/hr-ft^. Determine the heat transfer coefficient, 

13. Determine the coolant velocity required to keep the temperature 
rise across the convective layer below 65 F for the fuel and coolant 
conditions of Problem 12. 

14. In a nuclear fuel rod the volumetric thermal source strength varies 
axially as q' " (Z) = q^' " cos (nZ/H^) . The energy generated in this 
rod is removed by a coolant channel which has a mass flow rate of 

m and an inlet temperature of T^^. Designate the fuel diameter as D^ 
and the cladding OD as D^. Develop an expression for the axial 
variation of the temperature at the inside surface of the cladding. 
Your expression should contain only T^^, qQ*"* 

h, k , and Z. 

c 

15. From the expression developed in Problem 14 determine the location of 
the maximum cladding temperature in the reactor. 

16. Determine the location of the maximum fuel temperature for the fuel 
rod - coolant channel system of Problem 14. 

17. A 2760 Mw(t) PWR contains 31,000 fuel rods each composed of 0.380 in. 
diameter UO^ fuel pellets in 0.440 in. OD tubes of Zircaloy - 4 cladding. 
The fuel rods are spaced 0.580 in. in a square array. The reactor core 
height is 150 in. The coolant inlet conditions are 2250 psia, 553 F 

and the average mass velocity is 2.7 x 10° Ibm/hr-ft^. The thermal 
source strength varies as q**' = q^**' cos (HZ/H^) and H = 165 in. 
Determine the maximum fuel temperature and the maximum cladding 
temperature in the core . 

18. Determine the % reduction in coolant flow that will lead to incipient 
surface boiling for the reactor described in Problem 17. 

19. Determine the % overpower that the reactor described in Problem 17 
can produce without net steam formation. What fraction of the core 
experiences surface boiling at this limiting case? 
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20. A 2450 Mw(t) PWR contains 37,000 fuel rods with an active fuel length 
of 153 in. Each fuel rod consists of 0.362 in. diameter UO^ pellets 
enclosed in 0.420 in. CD Zircaloy - 4 cladding. The fuel rods are 
spaced 0«558 in. in a sqtiare array. The inlet coolant conditions are 
455 F, 2500 psia and the average velocity in the core is 15.8 ft/sec. 
The axial power distribution is a truncated cosine function. Assume 
an extrapolation length of 8 in. Determine the location and magnitude 
of the maximum fuel, cladding, and coolant temperatures. 

21. Determine the % of overpower that would result in incipient surface 
boiling in the core described in Problem 20. At what location does 
this occur? 

22. One of the coolant channels in the core described in Problem 20 is 
obstructed such that there is a 30% decrease in coolant flow. When 
and where does surface boiling occur? Does net void formation occur? 
If so, where? 
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